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Abstract. We present the A'-band (8-13 yum) spectrum of the hidden compact radio super-star cluster in NGC5253, C2, 
obtained with TIMMI2 on the ESO 3.6 m telescope. The spectrum is characterised by a rising continuum due to warm dust, a 
silicate absorption and a strong [S iv] line at 10.5 yum. Weaker lines of [Arm] at 9.0 yum and [Ne ii] at 12.8 /jm are also present. 
The continuum can be modeled by an optically thick emission from hot (T^ = 253 ± 1 K) dust emission extinguished by a 
cold foreground dust screen and a silicate absorption feature with A^-,] = 0.73 ± 0.05 mag. We show how the spatial scale of 
the observations greatly determine the mid-IR appearance of NGC 5253 and the important implications that this has on the 
interpretation of line fluxes in terms of the properties (age, IMF, etc.) of the embedded cluster. 

We have modeled the observed line fluxes towards C2 using photoionisation models with the most recent spectral energy 
distributions available to describe the integrated properties of the stellar cluster. The detailed dependence of the mid-IR lines 
on parameters such as the cluster age, upper mass cutoff and power law index of the IMF, as well as the local abundance, the 
presence of internal dust and the density structure is largely discussed. Strong constraints on the geometry based on high spatial 
resolution observations at different wavelengths - near-IR (HST and Keck), mid-IR (TIMMI2) and radio (VLA) - allows us to 
restrain the ionisation parameters to values log U > -0.5 dex. This constraint on U lead to two possible solutions for the age and 
upper mass cutoff of C2: 1) a young (< 4 Myr) cluster with a "non-standard" IMF with a low upper mass cutoff M„p< 50 Mq, 
and 2) a cluster of ~ 5 - 6 Myr with a standard high upper mass cutoff (Mi,p~ 100 Mq). A young age (< 4 Myr) would agree 
with the lack of supernovae signatures in C2 and in case of being confirmed, would be the first indication for a "non-standard", 
low upper mass cutoff of the IMF for an individual massive cluster. An older age of ~ 5-6 Myr would imply that it is possible to 
"contain" and hide such a compact cluster for a longer time that what it is generally thought. Arguments in favour and against 
these two scenarios are presented. 

The origin of the [O iv] 25.9 yum emission measured by ISO and the possible presence of an intermediate mass black hole inside 
C2 are also addressed. 

Key words. ISM: lines and bands - ISM: Hii region - Galaxies: starburst - Galaxies: star clusters - Galaxies: individual: 
NGC 5253 - Infrared: galaxies 



1. Introduction 

Near to mid-IR observations suffering from less or little ex- 
tinction are thought to provide more accurate spectral di- 
agnostics to quantify massive star formation and to distin- 
guish between stellar (s tarburst) and other (AGN) activity (cf. 
iGenzel & CesarskvlbOO O). With the near routinely access to 
mid-IR observations and thanks to their increased sensitivity, 
diagnostics in this spectral domain have been measured for a 
large number of Galactic and extra-galactic objects with the 
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Infrared Space Observatory (ISO). Furthermore, observations 
with the Spitzer Space Telescope have just begun. 

Despite the important advantage of strongly reduced ex- 
tinction, mid-IR satellite observations are potentially hampered 
by limitations due to their large aperture, especially for extra- 
galactic targets. As high spatial resolution mid-IR observations 
show, starburst galaxies and related objects can exhibit com- 
plex structures, numerou s "knots" etc. both in the continuum 
as in emission lines (e.g. iBoker et alJIlQQStlSoifer et alJl200ll 
l2002h on scales much smaller than the spectroscopic apertures 
(e.g. typically 280-380 arcsec^ for ISO/SWS). Furthermore, 
different diagnostics (e.g. PAH and continuum emission, high 
and low excitation fine structure lines) may originate from dif- 
ferent spatial regions. In consequence, it is by no means clear if 
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and to what extent such spatially integrated or "global" spectra 
can be used for various diagnostics purposes. 

For example, the reason for a relatively low average excita- 
tion in starburst galaxies as measured by the fine structure line 
ratio [Ne iii]/[Ne ii] (cf. Thornlev et al. 20o3) is not completely 
clear and has led to somewhat confl icting or inconsistent in- 
terpretations. IThornlev et aP tOOdt) argue e.g. that the low 
[Ne iii]/[Ne ii] ratios are due to short burst timescales and aging 
effects allowing, however, for a "normal" Initial Mass Function 
(IMF) including massive stars up to ^ 50-100 Mq. In contrast, 
[Rigbv & Rieke (2004) argue for a deficiency of stars with M ^ 
40Mq or suggest that massive stars or star clusters in massive 
galaxies spend a significant fraction of their lifetime hidden 
even at near to mid-IR wavelengths. Martin-Hernande z et alJ 
(l2002bl) . on the other hand, have pointed out a dependence of 
the observed [Ne iii]/[Ne n] ratio of starbursts on metallicity and 
argue that a proper treatment of the metallicity effects on stellar 
tracks, atmospheres and nebular models must be taken into ac- 
count in the comparisons. Finally one may wonder how reliable 
nebular models assuming a single "average" ionisation param- 
eter are and to which extent the low [Ne m]/[Ne ii] ratio is due 
to a non-negligible c ontribution of [N e ii] from a diffuse low ex- 
citatio n medium (cf . ISchaereJ2003l) . Although lThornlev et alJ 
(|20P0') estimate a diffuse contribution to be insignificant, this 
remains to be demonstrated observationally. 

In short, given the complexity of modelling and interpreting 
large aperture fine structure line emission and possibly other 
mid-IR diagnostics, there appears to be an important need for 
clarification based on simple objects representing, ideally, a 
single age and metallicity population with a well constrained 
geometry. The latter, in particular, is of prime importance to 
constrain the ionisation parameter of the nebular model and has 
generally not been treated with enough care. 

With these aims in mind we have recently started gath- 
ering ground-based high spatial resolution mid-IR observa- 
tions (imaging and spectroscopy) of several nearby star- 
bursts. Among them is the well studied bright starburst galaxy 
NGC 5253, a nearby Magellanic type irregular, observed at 
basically all wavelengths from X-ray to radio. Although this 
galaxy has a lower metallicity and high er excitation than the 
typical starbursts o bserved with ISO (cf. IThornlev et alJboOflt 
I Ver m a et al. 2003), the observations and analysis present here 
can provide interesting insight on the difficulties related to low 
spatial resolution spectroscopic observations and on the prob- 
lems discussed above. 

The galaxy NGC 5253 is also of prime interest for stud- 
ies related to super star clusters (hereafter SSCs), as it is 
one of the first starbursts where a very compact and hidden 
SSC was found from high spatial resolution radio observations 
jTurner et alJll998l) . Such optically thick bremsstrahlung ra- 
dio sources represent most likely the earliest phases of SSCs 
obser ved so far (e.g. iKobulnickv & Johnson 1999: Vacca et all 
I2OO2I) . As such, they provide invaluable insight on the prop- 
erties of SSCs, their formation and evolution. Furthermore, 
studying the youngest star forming regions yields in principle 
the most accurate constraints on the true upper limit of the IMF, 
since least affected by effects such as mass segregation, evapo- 
ration, and by the short finite lifetime of the most massive stars. 



Interestingly, the compact radio SSC in NGC 5253, often 
also referred to as the radio "supernebula", also stands out by 
its large contribution to the total int egrated flux of the g alaxy 
at various wavelength. For instance, 'Goriian et ajj ll200lh esti- 
mate that at least 15% o f the bolo metric luminosity of NGC 
5253 is due to this SSC. IVanzi & Sauvagd ( 120041) find that it 
dominates the total galaxy emission at wavelength longer than 
2-3 fj.m. Recently, HST/NICMOS observations have permit- 
ted resolving the region close to the supernebula revealing the 
presence of a n uclear double star cluster separated by 0'.'3-0'.'4 
( Alonso-Herrer oet alJl2004l) . The eastern near-IR cluster (CI) 
is identified with a young optical cluster, whereas the western 
star cluster (C2) coincides wi th the radio superne bula. Keck 
observations of NGC 5253 C2 dCioriian etalJ200ll) at 11.7 fim 
reveal an unresolved source with a diameter of 0'.'58 and a total 
flux of 2.2 Jy. R adio observations at mm and cm wavelengths 
( Tur ner et alJ20 00: Turner & Beck 2004) show a slightly elon- 
gated source of about ~ Of! 1 xO'.'05. As we will see below, these 
high spatial resolution observations combined with our mid- 
IR spectroscopy provide, by means of detailed photoionisation 
models, a unique opportunity to study in detail the properties 
of this extreme hidden massive star cluster 

The core of the paper is structured as follows. Our observa- 
tions are described in Sect.|2| Immediate results from our spec- 
tra are discussed in Sect. |3| Detailed photoionisation models 
for the compact SSC C2 are presented in Sect.|3 Our principal 
results on the age, initial mass function and other properties of 
the SSC, as well as various implications are discussed in Sect. 
|5| The main conclusions of the paper are summarised in Sect. 

m 

2. Observations and data reduction 

Our new infrared data on NGC 5253 were obtained as part of a 
programme with the Thermal Infrared MultiMode Instrument 
(TIMMI2) on the ESO 3.6 m telescope (La Silla Observatory, 
Chile) to observe young starburst galaxies. 

The A^-band spectrum of NGC 5253 C2 was obtained on 
2003 March 2 1 . We used the 1 fim low-resolution grism which 
ranges from 7.5 to 13.9 jum and has a spectral resolving power 
A/AA ~ 160. The slit width used was l'.'2x70", with a pixel 
scale of 0'.'45. The slit was positioned across the bright infrared 
supernebula C2 and in the north-south direction. In order to 
correct for background emission from the sky, the observations 
were performed using a standard chopping/nodding technique 
along the slit (where the object is observed at two different po- 
sitions on the slit) with an amplitude of 10" in the north-south 
direction. The effective exposure time (on-source) was 32 min. 
The observation was performed at an airmass of ~ 1.2. The star 
HD 123139 was observed right before and after C2 and served 
as both telluric and flux standard star This is a primary ISO cal- 
ibration standard star and is described in detail at the TIMMI2 
webpage'. The synthetic calibrated spectrum for this standard 
star is given bv lCohen et alJ ( Il999l) . 

The data processing included the removal of bad frames 
and the coaddition of all chopping and nodding pairs. This left 

' www.ls.eso.org/lasilla/sciops/timmi 
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Fig. 1. (a) Top panel shows the A^-band spectrum of the embedded cluster NGC 5253 C2. The spectrum is characterised by a 
rising dust continuum and by the atomic fine-structure lines of [Ar iii] at 9.0 //m, [S iv] at 10.5 jim and [Ne ii] at 12.8 //m. Bottom 
panel shows the statistical uncertainty at each wavelength, (b) Top panel shows the best optically thin and optically thick models 
for the dust continuum emission (see Sect. 13. 3> . Bottom panel shows the continuum- subtracted spectrum when applying the 
optically thick model for the continuum. 



us with one single image with one positive and two negative 
long-slit spectra. These were combined with a simple shift- 
and-add procedure which increased significantly the signal-to- 
noise. The spectra of the target and standard star were ext racted 
using the optimal extraction procedure developed by Hornd 
( Il986<) . ideal for unresolved point sources. This procedure ap- 
pUes non-uniform pixel weights in the extraction sum in or- 
der to reduce the statistical noise in the extracted spectrum to a 
minimum while preserving its photometric accuracy. Slit losses 
are negligible, since the full-width-at-half-maximum (FWHM) 
along the spatial direction for both target and standard star is 
r.'l, smaller than the slit width (r.'2). The calibration of the 
spectroscopic data included 7 ) the removal of the telluric fea- 
tures, which was done by dividing by the spectrum of the stan- 
dard star, 2) the removal of the spectral features of the standard 
star and 3) the absolute flux calibration. These last two steps 
were conjointly done by multiplying by the synthetic spectrum 
of the standard star. Wavelength calibration is straight forward 
since a table with the pixel-to-wavelength correspondence is 
provided by the TIMMI2 webpage. 



3. Results 

3.1. The N-band spectrum 

Figure ^a shows the spectrum of the embedded cluster 
NGC 5253 C2. Only the vaUd range between 8.15 and 13 i^m 
is shown. The spectrum is characterised by a rising continuum 
due to warm dust, a silicate absorption around 9.7 yum and the 
fine-structure fines of [Arm] at 9.0 //m, [S iv] at 10.5 fim and 



[Ne ii] at 12.8 yum. Tfie [Ne ii] line is very weak. The presence of 
silicate absorpt ion contradicts the best fitting model of the nu- 
clear SED bv .Vanzi & Sauvag3(l2004l) . which does not require 
a sificate component. T he ISO/SWS spectrum of NGC 5253 by 
ICrowther et alJ ( 1 19991) shows no evidence of absorption at the 
location of the silicate band. However, the silicate absorption 
might be hampered by the large ISO aperture (~ 14"x20") and 
by the low signal-to-noise of the spectrum. A direct compari- 
son with the ISO/SWS spectrum also shows that the continuum 
flux level of NGC 5253 C2 is about 3 times lower that that ob- 
served by the large ISO aperture. 

The spectrum does not display PAH emission features. 
This is confirmed by the L-band spectrum of the nuclear re - 
gion of NGC 5253 obtained by .Alonso-Herrero et alJ ( |2004 . 
which shows the non-presence of the 3.3 fim P AH feature. The 
ISO spectrum integrated over the whole galaxy dCrowther et alJ 
1222) shows no signature of PAH emission features neither. 
Most massive star-forming galaxies show emission bands at 
3.3, 6.2, 7.7, 8.6 and 11.2 /im, generally attributed to emis- 
sion by Polycyclic Aromatic Hydrocarbon (PAH) molecules 
(e.g. lAHamandola et all Il989l) . Hence, the non-detection of 
PAH features towards NGC 5253 C2 mak es it look more like 
an AGN than a star forming reg ion (e.g. iGenzel & Cesarskvl 
l2000t ISiebenmorgen et al]l2004i . However, it is shown from 
ISO observations that low metallicity dwarf galaxies such as 
NGC 5253 tend to show weak of no mid-IR P AH emission fea- 
tures (e.g. iThuanet al.ll 19991 iMaddenlEool . Moreover, PAH 
features are not present in the mid-IR s pectra of highly ern - 
bedded massive star formation sites (e.g. IPeeters et al.ll2004l) . 
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The effects of low metallicity and a harsh environment in 
NGC 5253 C2 probably suppress the formation of PAH emis- 
sion features. 

3.2. Line fluxes observed at different spatial 
resolutions 

Line fluxes in the A^-band spectrum were measured by fitting a 
Gaussian and are Usted in Tabled We have compared these line 
fluxes with IR measurements for NGC 5253 from the literature. 

lRecketai1lll994 observed the nuclear region of the galaxy 
with a resolution of 1'.'6 and obtained fluxes for the [Arm], 
[S IV] and [Neii] lines of, respectively, 3.6 x lO^^o^ 39 ^ ^q-^o 
and < 1 X 10"^'^ Wcm"^. The [Arm] line flux is comparable to 
the one we have measured with TIMMI2, while the line flux of 
[S iv] is about 2 times higher The reason of this discrepancy is 
not clear to us. 

The central region of NGC 5253 has also been observed by 
the SWS on b oard ISO, resulting in a spectrum ranging from 
2.3 to 45 yum JCrowtheret alJI 19991) . The ISO/SWS apertures 
are approximately 14" x 20" up to 12 fim, and 14" x 27" be- 
tween 12 and 19.6 fim. The line fluxes measured by ISO/SWS 
are (5.9 + 0.5) x 10"^°, (34 + 3) x 10"^" and (8 + 1) x 10"^° 
Wcm^^, respectively for the [Arm], [Siv] and [Neii] lines 
JVerma et alJl2003l see also Tabled, where we have adopted 
as uncertainties the typical ISO/SWS Icr absolute flux accu- 
racy (Peeters et al. 2002). Line ratios between the ISO/SWS 
and TIMMI2 fluxes are 1.3 + 0.3, 1.9 + 0.1 and 5.7 ± 2.0. Ar++, 
S^^ and Ne^ have ionisation potentials of 27.6, 34.8 and 21.6 
eV, respectively. Hence, in view of the comparison between 
the line fluxes measured towards NGC 5253 C2 by TIMMI2 
and those measured by the large ISO aperture, we see that a 
large fraction of the emission from the high-excitation lines 
([Arm] and [S iv]) in the central starburst region of NGC 5253 
comes from the compact and embedded H 11 region associated 
to NGC 5253 C2. However, only -20% of the [Neii] emission 
comes from this H 11 region. Other sources such as the remain- 
ing observed clusters and the diffuse ISM must contribute to 
the [Ne 11] line flux measured by ISO/SWS. 

Ratios of fine-structure lines are commonly used in a di- 
versity of diagnostics. For instance, the ratio of the high- 
to-low-excitation emission lines [Oiv] 25.9/[Neii] 12.8 fim is 
used, in combination with the 7.7 fim PAH line-to-continuum 
ratio, to separate AGN dominated and star formation domi- 
nated nuclei (cf. Genzel et al. 1998)). As well, ratios such as 
[Nem] 15.5/[ Nenl 12.8 //m are used to infer the age of the 
starburst (e.g lThornlev et al1l2000 ). The presence of IR-bright 
hidden clusters such as the one observed in NGC 5253 and 
their contribution to the global atomic line emission might thus 
obfuscate diagnostic diagrams used to determine the ultimate 
physical process powering galactic nuclei and studies of the 
IMF. 

We also compare the Bra line fluxes measured by 
iTurner et all (2003) using NIRSPEC on the Keck Telescope- 
with an aperture of 0.579" centred on NGC 5253 C2 - and by 
ISO/SWS. These fluxes are, respectively, (3.4+ 1 .0)x 10"^" (see 
Table IB and (4.5 + 0.4) x lO^^o Wcm^^ (.Verma et al...2003.) . 



Both line fluxes are in good agreement, indicating that practi- 
cally all the Bra emission is produced by NGC 5253 C2. 

3.3. Extinction 

3.3.1. Extinction in tine near-IR 

iTurner et alJ ( l2003l) observed the Brackett y (2.17 //m) 
and a (4.05 fim) recombination lines of hydrogen towards 
NGC 5253 C2 using NIRSPEC on the Keck Telescope. The ob- 
served and extinction-corrected line fluxes are listed in Tabled 
The lines were corrected for extinction considering that the 
intrinsic Fa/Fy flux ratio is 2.8 (for an^electron temperature 
Te = 12 000 K; Hummer & Storevlll987ll) and an extinction 
law Ai oc A^^''' (Mathis 1990). The value of the electron tem- 
perature is deter mined for visible Hii regions in NGC 5 25 3 
(Ca mnbell et all 1 19861: IWafsh fc Rovlll989t iKobulnickv etalJ 
11997). We have recomputed the extinction values obtained by 
[Turner et al. (2003) in order to include a detailed error analy- 
sis. The extinctions with their respective uncertainties are: Ay - 
1.9+0.6 mag, A„ = 0.7+0.2 mag and AK(2.2Aim) = 1.84+0.63 

mag. 

Assuming that Av/Ak ~ 9.3 forRy = 3.1 ( lMathislll99(il) . 
we obtain that NGC 5253 C2 is extinguished by 17 + 6 mag in 
the optical V-band. This value of the extinction agrees w ell with 
the near-IR estimate bv lAlonso-Herrero et alJ ll2004h . 18+5 
mag in an aperture of 0'.'6x3". The modelling of the near- 
infrared to millimetre SEP by l^nzi & Sauvage (2004) and 
lAlonso-Herrero et alJ ( l2004h gives comparable values for the 
extinction (16-21 mag). We note that this extinction range is 
obtained by Vanzi & Sauvage (2004) in their alternate model 
which does not consider the HST optical data. Dropping the 
optical flux as a constraint on the dust modelling is more ap- 
propriate as the dominant IR emitting region (C2) and the opti- 
cally emitting regions are unrelated. In contrast, their principal 
model, which wrongly includes optical data, gives a lower ex- 
tinction of ~ 8 mag. 

3.3.2. Extinction in tine A'^-band 

In order to determine the intrinsic line fluxes in the A^-band, we 
need to correct for extinction, which, in the mid-IR, is mainly 
caused by silicate absorption. Following Okamoto et al. (200l[ 
l2003t) . we have adopted a simple model to derive the extinction 
in the A^-band, where the emission from hot dust in or at the 
outer edge of the ionised region is assumed to be extinguished 
by cold foreground dust along the line-of-sight. This represen- 
tation of the dust distribution is supported by the UV study of 
Meurer et alJ lfT995). who concluded that the dust geometry to- 
wards NGC 5253 could be described in terms of a foreground 
dust screen near the starburst. In such model, the emission of 
the dust at a wavelength A is described as: 

F, = AB,(rd)e,e-^' , (1) 

where A is a scaling factor, is the temperature of the emitting 
hot dust, Q is the dust emissivity and t,i is the optical depth. We 
assume that the foreground dust has the same grain properties 
as the emitting hot dust and hence, ta can be written as TsueA, 
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[Arm] 9.0 


4.4 ± 1.0 


7.6+1.8 


0.60 +0.04 


1.2 


[Siv] 10.5 


18.0 + 0.7 


29.9+1.6 


0.55 +0.04 


1.2 


[Neil] 12.8 


1.4 + 0.5 


1.7 + 0.6 


n T2 _i_n m 
U.ZJ +U.UZ 


1 

l.Z 


Turner et al. (20031: 








Bry 2.17 


0.4 ±0.1 


2.2+1.4 


1.9 + 0.6 


0.579 


Bro- 4.05 


3.4+1.0 


6.2 + 2.2 


0.7 + 0.2 


0.579 


Lumsdenetal. (1994): 








He I 2.06 


0.240 ± 0.004 


1.6+ 1.0 


2.1+0.7 
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Selected mid-IR line fluxes measured bv ISO/SWS fVerma et al. 2003): 




[Am] 7.0 


2.3 + 0.2 


2.6 + 0.2 


0.13+0.01 


14x20 


[Arm] 9.0 


5.9 + 0.5 


10.2 + 0.9 


0.60+0.04 


14x20 


[Siv] 10.5 


34 + 3 


56 + 5 


0.55+0.04 


14x20 


[Neil] 12.8 


8 + 1 


10+ 1 


0.23+0.02 


14x27 


[Neiii] 15.5 


28.9 + 4.6 


34 + 5 


0.19+0.01 


14x27 


[Siii] 18.7 


15.9 + 2.5 


20 + 3 


0.25+0.02 


14x27 


[Oiv] 25.9 


0.9 + 0.1 


1.0 + 0.1 


0.14+0.01 


14x27 



" Before extinction correction. * After extinction correction. 



where Tsu is the peak optical depth of the silicate. We have 
considered two extreme cases for the optical emission of the 
emitting hot dust: an optically thick approximation (where ca = 
1) and an optically thin approximation (where oc t,0. The 
adopted is based on a mixture of carbonaceous grains (PAH 
and graphite- like) and amorphous silicate g rains with Ry - 3.1 
developed bv lWeingartner & Drainel(l200ll) which successfully 
reproduces the observed Galactic interstellar extinction. 

The 8-13 jum spectrum was therefore fitted using the above 
model, which has 3 free parameters: A, and t^^. We per- 
formed least-squares fits to the data. The best fits for the op- 
tically thick and optically thin approximations are shown in 
Fig. n]b. We obtained a significantly better fit using the opti- 
cally thick approximation (with a resulting rtductd-x^ equal 
to 1.5) than using the optically thin approximation (reduced- 

- 4.8) and therefore, we will only consider the first case 
hereafter Figuref^b shows the continuum-subtracted spectrum 
using the optically thick model. 

The best fit (optically thick approximation) gives A = (4.2+ 
0.1) X 10"'^ (with Ba given in Wcm-Vm"'), = 253 + 1 K 
and Tsii - 0.67 + 0.05. The error analysis was done by follow- 
ing how changes as the parameters in the model vary. The 
final errors in the best-fit parameters were calculated by deter- 
mining how much the parameters have to change to get to the 
1% probability level. 

The values of the extinction estimated at the wavelengths 
of the A^-band fine-structure lines are given in Table ^ to- 
gether with the line fluxes after corrected for extinction. The 
extinction in magnitudes is related to the optical depth through 



Aa — 1.086t^. We also give the extinction estimates and cor- 
rected line fluxes of the mid-IR lines measured by ISO. These 
line fluxes must be taken as upper limits because of the large 
ISO aperture. 

We can now draw the complete 2-13 jum extinction law 
for the embedded cluster, where the range between 2 and 
-7.5 jum is given by Aa = (1.84 + 0.63)(A/2.2 fim)-^ '^ (see 
Sect. 13 3. and the range beyond is given by the silicate pro- 
file bv lWeingartner & Draind ( l200ll) wifli A^u = 0.73 + 0.05 . 
Since we have obtained that Ay = 17 + 6 mag (see Sect. 13. 3.11 . 
we have that Av/Asn = 23 + 8, in good agreement with the 
Av/Asii ratio of 18.5 ± 1.5 found bv lRoche & AitkerJ ( ll984 for 
the local diffuse interstellar medium. 

3.4. Ionic abundances 

Ionic abundances can be determined from the measured 
strengths of the fine-structure and H i recombination lines. The 
ionic abundance of a certain ion X^' with respect to hydro- 
gen (X^' /H^) can be de termined using the following expression 
(e.g.lRubinetalJl988l): 

where Fx*i and Fm are the extinction-corrected fluxes of any 
line produced by X^' and H i, and ex+i and ehi are their respec- 
tive emission coefficients. This expression assumes that 1) the 
nebula is homogeneous with constant electron temperature and 
density, 2) all the line photons emitted in the nebula escape 
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without absorption and therefore without causing further up- 
ward transitions and 3) the volume occupied by X^' and is 
the same. 

The emission coefficients depend on Tp, Wp a nd the relevant 
atomic parameters ( Martin-Hernande z et al ]|B02a). Adopting 
re=12 000 K (see Sect, imi and «e=5 x lO'* cm^^ (see 
Sect. 14. 1> . we obtain qAriii] 9.0 = 8.80 x 10"^', qsiv] 10.5 = 
1.41 X 10-20, e[Neii] 12.8 = 7.27 x 10-^2 and eBr« = 7.60 x 
10^2^ erg cm-* s"'. Using these values for the emission coef- 
ficients and the line fluxes listed in Table \l\ the ionic abun- 
dances with respect to we estimate for Ar^^, S^^ and Ne^ 
are, respectively, (1.0 + 0.4) x 10"^ (2.6 ± 0.9) x lO"*" and 
(2.9+ 1.4) X 10-^ 

Finally, we can estimate upper limits for the elemental 
abundances of Ar, S and Ne using the [Am] 7.0, [S iii] 18.7 and 
[Neiii] 15.5 fim line fluxes measured by the larger ISO aper- 
ture. The emission coefficients of these lines for re=12000 K 
and «e=5 x 10"* cm"^ are e[Arii] 7.0 = 1-22 x 10"^", e[siii] 18.7 = 
2.09 X 10^2' and qNeiii] 15.5 = 1-24 X 10"^' erg cm^ s'K Hence, 
upper Umits to the ionic abundances with respect to H^ of Ar^, 
S^^ and Ne++ are, respectively, 3.5 x 10"^, 1.6 x 10"^ and 
4.7 X 10"^. Taking into account these values, we can give upper 
limits to the total elemental abundances of Ar, S and Ne assum- 
ing that most of the argon and neon is singly and doubly ionised 
and that most of the sulphur appears as S^^ and S^^. The up- 
per limits we obtain for Ar/H, S/H and Ne/H are, therefore, 
1.7 X 10"*, 1.9 X 10"^ and 5.1 x 10"^, which are approximately 
0.7, 0.9 and 0.4 times the solar elemental abundances. The up- 
per limit that is likely closer to the true elemental abundance is 
that of Ne/H since most of the [Ne iii] line flux must come from 
the H II region formed by the embedded star cluster C2 (as it 
is the case of the [S iv] fine flux, cf. Sect. 13. 2> . The upper limit 
obtained for Ne/H (O.4x[Ne/H]0) agrees wefl with the global 
metallicity of the galaxy (Z ~ 0.302©, see Sect. 14. U . 

4. Modelling the nebular lines 

We begin with a description of the observational constraints 
and input parameters for the modelling of the nebular lines. 

4.1. Observational constraints 

There are several observational facts that can be used to con- 
strain the photoionisation models: 



cally constant during the first 6 Myr after the burst of star 
fo rmation, with a value of th e order of the EW observed 
bv lAlonso-Herrero et alJ ('2004). decreasing drastically beyond 
this age. Hence, based on the observed EW, we can only ob- 
tain an upper limit for the age of the cluster of 6 Myr. This 
upper limit is independent of the upper mass cutoff" of the clus- 
ter initial mass function since variations due to different upper 
mass cutoffs are smaller than the quoted observational error. 
Indication of an early age comes as well from the fact that 
the nuclear radio and millimetre continuu m emission are al- 
most entirely due to thermal emission re.g. lBecketal]ll996t 
iTurner et a l. 2000), which rules out supernova remnants or ra- 
dio supernovae as the source of the emission and place an up- 
per limit of about 4 Myr for the age of the cluster Finally, 
the finding of optica lly thick thermal radio emission from the 
C2 region (iTurner et al...l998i) interpreted as an "ultra-dense" 
H n region tends to indicate even younger ages (^ 1 Myr; 
cf. lKobulnickv & Johnso n 1999:' V acca et al.^ 2002') as deduced 
from pressure/lifetime arguments. However, this argument can 
possibly be relaxed if the cluster mass is large enough al- 
lowing for its own gravitational confinement, as suggested by 
iTurneret alJ ( l2003l) . 

Geometry Recent 7 mm co ntinuum observations made with 
the VLA including Pie Town jTurner & Beck 2004) shows that 
the supernebula has a size of 99 + 9 mas by 39 + 4, with a 
position angle of 6 ± 4deg, or 1.6 pc by 0.6 pc for the above 
distance of 3.25 Mpc. The nebula is smaller than most extra- 
galactic super star clusters (cf. Table|5}. 

Chemical composition Optical and ultraviolet spectroscopy 
obtained with the HST Faint Object Spectrograph at 
three locations in the central Hii complex of NGC 5253 
(Ocobulnicky et al. 1997) gives an oxygen abundance of 12 -1- 
log(0/H) - 8.16 + 0.06 i n agreement \yith pr evious ground- 
based spectroscopy (e.g. IWalsh & Rovl Il989l) . Assuming a 
solar oxygen abundance of 12 H- log(0/H)Q - 8.69 
jAUende Prieto et alJl200ll) . NGC 5253 is a metal-poor galaxy 
with a metaflicity Z ~ 0.30 x Z© = 0.006. This metaflicity 
agrees well with the upper limits for the elemental abundances 
of NGC 5253 C2 derived in Sect. 13.41 in particular, it is close 
to the upper limit obtained for Ne/H, which is close to the true 
Ne abundance of the embedded cluster (see also Sect. l3.4t . 



Distance The dista nce to NGC 5253 is w ell known from its 
Cepheids. Recently. iFreedman et al.l ( I2OOII) has derived a dis- 
tance modulus of fiQ = 27.56 + 0.14 from 17 Cepheid can- 
didates. This distance modulus is equivalent to a distance of 
3.25+0.21 Mpc. 

Age Recently, lAlonso-Herrero et 

Z] 1I2OO4I1 have estimated 
the age of the stellar population of NGC 5253 C2 based on the 
equivalent width (EW) of Pao- (log(EW(Paa))=3.48 + 0.23). 
They derive an age of 3.3 ± 1.0 Myr. We have, however, com- 
puted the EW of Pao- predicted by our photoionisation mod- 
els (Sect. 14.21 and thereafter). The predicted EW is practi- 



Electron density The electron density can be constrained from 
radio continuum observations. Subarsecond resolution obser- 
vations of NGC 5253 C2 indicates that t he region is par- 
tially optically th ick at cm wavelengths jBecket alJll996l: 
iTurner et alJl200ol) . However, it has been recently observed at 
7 mm (43 GHz) bv lTurner & Beck (2004). The nebular emis- 
sion is optically thin at this wavelength. The nebula contains 
8.5 mJy of flux at 7 mm. Assuming that the line-of-sight dimen- 
sion is equal to the transversal size of the region (~1 pc if we 
take the geometrical average of the major and minor axes of the 
core as a representative size, see above), the rms electron den- 
sity (calculated using Eq. A. 2.5 bv .Panagia & Walmsley..l978< 
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) is ne~ 5 X 10"* cm"^. This electron density is typical of ul- 
tracompact H ii regions and it is high for an H ii region of this 
size. Densities and sizes for Galactic H ii regions follow a rather 
tight relation of the form log(«p) ~ 3 - 0.9 x log D, wher e D is 
the di ameter in pc dOarav et alii 19931 iMar tm-Hernandez et alJ 



Table 2. Observational constraints. 



120031) . A Galactic Hii region of about 1 pc in diameter would 
have an electron density of the order of 10^ cm"^. 

H I recombination lines Recombination lines of hydrogen de- 
pend on the rate of Lyman ionising photons, Qq. The observed 
and extinction-corrected line fluxes for the Brackett y and a 
lines (see Sect. 13. 3.11 are listed in Tabled 



He I recombination lines 'Lums den et alJ (1 1994*) observed the 
central 3"x3"region of NGC 5253 in the K-band. They mea- 
sured a Bry flux equal to 0.5 x 10"^*^ W cm~^. The good a gree- 
ment with the Bry flux measured by Turner et"al] 12003') (see 
Table Q indicates that the i onised emission in the K-band is 
dominated by NGC 5253 C2. lLumsden et alJ ( 1 1994 also mea- 
sured the 2.06 jum (2'P - 2^5) Hei line in this band. The ob- 
served and extinction-corrected line fluxes of this He i line are 
listed in Table The ratio He 1 2.06/Br7 depends on the ratio 
of the He (hv > 24.6 eV) to H (hv > 13.6 eV) ionising photons. 

Mid-infrared fine-structure lines Mid-IR fine-structure 
lines can be used to probe the stellar SED in the EUV 
range. The observed and extinction-corrected line fluxes of 
the mid-infrared lines observed by TIMMI2 ([Arm] 9.0, 
[Siv] 10.4 and [Neii] 12.8 fim) are listed in Table [T] We 
also list other mid-infrared line fluxes measured by ISO 
([Am] 7.0, [Neiii] 15.5, [S iii] 18.7 and [Oiv] 25.9 jum), which 
must be taken as upper limits because of the large ISO aperture. 

A summary of the above observational constraints is shown 
in Table|2l 

4.2. Input parameters and assumptions 

We compute sets of nebular models with the photoionisation 
code CLOUDY^ (Ferland et al. 1998) version 96.00-beta 4 us- 
ing MICE^ die IDL interface for CLOUDY created by H. 
Spoon. The computation is performed for a static, spherically 
symmetric, ionisation bounded gas distribution with an inner 
cavity. We assume that the gas is uniformly distributed in small 
clumps of constant density over the nebular volume and occu- 
pies a fraction e of the total volume. The shape of the stellar ra- 
diation field, the ionising photon luminosity (Qo), the electron 
density {rig.), the inner radius of the shell of ionised gas (Rin), 
the filling factor (e), and the chemical composition (Z) are the 
input parameters for the photoionisation models to predict the 
intensity of the emission lines. 

We have used the evolutio nary synthesis code 
Starburst99'* dLeitherer et al ] ll999ll version 4.0 to model 



see http : //thunder. pa. uky. edu/cloudy/j 



2 

^ see http://www.astro.rug.nl/~spoon/mice.html 
see ,http://www.stsci.edu/science/starburst99/. 



Parameter 


Value 


Bra 


(6.2 ± 2.2) X 10"-" W cm-2 


Hei/Bro- 


0.3 + 0.2 


[Ariij/Bro- 


< 0.6 


[Ar iii]/B ra- 


1.2 ±0.5 


re ivl/Rrn' 


4.8 + 1.7 


[Neii]/Bra 


0.3 ±0.1 


[Neiii]/Bra 


< 7.7 


[S mj/Bra 


<4.5 


^out 


0.8 pc 




5 x 10" cm-3 


z 


0.006 



the integrated properties of the stellar cluster This code is 
based on stellar evolution models of the Geneva group and 
uses enhanced mass loss tracks for masses above 12 Mo 



(Mevnet et al. 


1994) and standard mass loss tracks between 


0.8 and 12M^(Schaller et alJ ll992HSchaerer et alJ Il993albl: 


iCharbonnel et al 


.i.l993i). Starburst99 follows the evolution in 



the H-R diagram of a stellar population whose composition is 
specified by a stellar initial mass function (IMF). At a given 
age, the integrated SED is obtained by summing over the con- 
tributions of all stars prese nt and is built using 1 ) t he non-LTE 
WM-Basic steUar models jPauldrach et alJl200lt ISmith et al] 
i2002i) for O stars, which take into account the effects of stellar 
winds and line blanketing, 2) the f ully line-blanketed mo dels 
calculated with the CMFGEN code dHiUier & Miller' 1998) for 
W olf-Rayet (WR) stars and 3 ) the plane-parallel LTE models 
by iKuruczl ( 1 19931) for the remaining stars that contribute to 
the continuu m. We assume an instantaneous burst of star 
formation, a ISalpeteJ ll 19551) IMF with exponent a - 2.35 
(dN/dlnm oc m'""), a lower mass cutoff Mio„ = 1 M0 and 
an upper mass cutoff /V/yp set to 30, 50 and 100 Mq. We also 
assume that the stars evolve from the main sequence following 
the Z = 0.004 and Z = 0.008 steUar tracks (stellar tracks with 
Z = 0.006 are not currently available). We present models 
every 0.5 Myr for 10 Myr after the burst of star formation. 

The chemical composition of the gas is set to Z = 0.30 x 
Zq - 0.006 (see Sect. 14. 1> . We have adjusted the helium 
abundance Y according to Y - Yp + (A Y/AZ)Z, where F p = 
0.24 is the primordial helium abunda nce (lAudouzd Il987h and 
(AF/AZ) = 3 is an observed constant (Paget 1992). To arrive at 
the appropriate metal abundances, we have simply scaled the 
solar values. 

A change of the rate of Lyman ionising photons, electron 
density, inner radius of the shell and/or filling factor is equiv- 
alent to a change of the ionisation parameter (U) defined by 
CLOUDY as: 



(3) 



where is the distance between the ionising cluster and the 
illuminated surface of the gas shell and riu is the total hydro- 
gen density (ionised, neutral and molecular). Basically, U is 
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the dimensionless ratio of hydrogen ionising photons to to- 
tal hydrogen density. We fix an electron density equal to the 
rms electron density derived from radio continuum observa- 
tions (5 X lO'* cm"^, see Sect. 14. H . Qq is fixed individually for 
every photoionisation model so that the Bm line flux is repro- 
duced within 5%. Test models give that the resulting Qq varies 
slightly with the cluster age, with values between 1.6 x 10^^ 
and 2 x 10^^ s"' , in fair agreement with the ~ 1 x 10^^ s~^ esti 



mated from the 7 mm observations presented bv lTurner & Beckl 

tooA 

Once we have fixed rig and Qo, there is still a degeneracy 
between the inner and outer ra dius of the H ii region and its 
filling factor. From lOsterbrockl ( fl989) . we have that: 



Qo 



47T 



(4) 



where the rate of ionising photons emitted by the cluster (Qo) 
just balances the total number of H i recombinations to excited 
levels within the ionised volume 4jt(RI^^ - )/3, where /?out is 
the outer radius of the gas shell. The parameter q-b is the case 
B Hi recombination coefficient to all levels > 2. Using Eq.|4] 
we can express U (Eq. as a function of R^^, Rom and e as 
follows: 



Q'B 

U = — Hef 

3c 



R^ 



3 \ 



(5) 



where we have assumed that riu - «e- The size of the source ob- 
served byEmir^Bi£3 lEooi) is 1.6 X 0.6 pc (see Sect.lO. 
Hence, Rout cannot be larger than 0.8 pc. The outer radius of the 
model comes defined by the position where the ratio of electron 
to total hydrogen densities reaches a value of 0.90, i.e. near 
the H^-H" ionisation front. Ri„ and e, however, cannot be con- 
strained from observations, but we can set them to reasonable 
values. 

Equation0]can be re-written as: 



(1_,3)^3 ,,^6^, 2.14X10^^ ^^ 

Amtaa nt nz 



(6) 



where 7?out is in pc and 7?i„ - xRom, with < x < 1 . We have 
adopted Qo - 2x 10^^ s '. Physical solutions of this equation 
are those that satisfy the condition e < 1, which is equivalent 
to 



Ji:< 1 



y6 



R^ 



1/3 



(7) 



This condition leads to x < 0.75 for «e=5 x 10"* cm"^ and Rom - 
0.8 pc. We note that this condition does not allow values of 
^out < 0.7 pc and «£< 3.7 x 10'* cm"^. The minimum value 
allowed for the filling factor is e(x — » 0) = 0.6. Hence, the final 
constraints on the geometry are R^at - 0.8 pc, R^^ < 0.6 pc 
and 0.6 < e < 1. We remind that these three parameters are 
not independent, but have to satisfy the condition expressed by 
Eq.|6l 

We evaluate now how the above variations of and e 
modify the ionisation parameter U. Fig. |2l shows the varia- 
tion of U for an Hii region with Qo = 2 x 10^^ s ' and 



3 

cn 
O 




Fig. 2. Variation of ionisation parameter U with inner radius 
(Rin) for an Hii region with Qo - 2 x 10^^ s"' and «e=5 x 
10"^ cm"^. Ri„ - 0.6 pc is the maximum value allowed for the 
inner radius of an H n region with Rom = 0.8 pc. 

Table 3. Input parameters of the fiducial model. 



Parameter 


Value 




0.6 pc 


^out 


0.8 pc 


e 


1 


"e 


5 X 10* cm-3 


7 


03Zq 



ne=5 X 10^ cm"-' (Eq. |3}- ^in = 0.6 pc (which coiTesponds 
to a filling factor e = 1) is the maximum value allowed for the 
inner radius if Rom = 0.8 pc. The ionisation parameter for a 
nebula with Ria = 0.6 pc and e = 1 is, according to the figure, 
\ogU - -0.5 dex. However, regions with smaller inner radii 
will have larger ionisation parameters. 

We will adopt a reference geometry with Ri„ - 0.6 pc and 
e = 1 for the main bulk of models presented hereafter. A sum- 
mary of the input parameters of these fiducial models is pre- 
sented in Table |3l The effect of varying U (i.e. Rm and e) and 
other parameters such as the slope of the IMF, the nebular metal 
content, the internal dust content and density structure will be 
discussed in the following section. 



4.3. Results from the photoionisation models 

Figures |3] and 0] show the main results of the photoionisa- 
tion models compared to the (extinction-corrected) observed 
values listed in Table [2 Figure |3] shows the variation of 
the seven independent constraints from Table [2 ([Ar n]/Bra, 
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[Ariii]/BrQ', [Siii]/Bra', [Sivl/Bra, [Neii]/BrQ', [Neiii]/Bra' and 
Hei2.06/Bra') as a function of time since the burst of star for- 
mation. We note that these models reproduce the observed 
Bra Hne flux within 5%. Figure 0] shows the variation of 
the mid-infrared hne ratios [Siv]/[Neii], [Ar iii]/[Ne ii] and 
[Neiii]/[Neii]. 

In agreement with ear lier combin ed starburst and 
photoionisation models (e.g. ISchaerer & Stasihska 1999a b; 
iRigbv& Riekl \2Q04), we observe that ratios such as 
[Siv]/[Neii], [Arm]/[Neii] and [Ne iii]/[Ne n] (see Fig.gl de- 
crease smoothly during approximately the first 3 Myr as the 
hottest O stars leave the main sequence, increasing again from 
3 to 4-6 Myr as a consequence of the contribution of the WR 
stars and the remaining main sequence stars to finally drop 
afterwards. This clearly reflects the importance of WR stars, 
which contribute to the ionisation of the nebula after the hottest 
O stars have left the main sequence. This behaviour is espe- 
cially noticeable in the case of models with Myp = 100 Mq. 
This contribution of the WR stars is also observable in Fig.|3] 

These figures clearly show the dependence of these line ra- 
tios on the stellar content and metallicity. A higher Mup implies 
a harder stellar spectrum and consequently, a higher degree of 
excitation. Similarly, a lower stellar metallicity leads as well to 
a harder stellar spectrum. 

We will focus first on the subsolar predictions. The compar- 
ison between the results of these photoionisation models and 
the (extinction-corrected) observed line ratios clearly shows 
that it is not possible to constrain at the same time both the 
age and Mup of the stellar cluster. However, some conclusions 
can be derived. First, it is shown that NGC 5253 C2 has an age 
of less than 6 Myr. All model predictions beyond this age fall 
far from the observations. A cluster with Mup = 100 Mq repro- 
duces well all the observed lines for an age of ~ 5 - 6 Myr. 
Second, a closer inspection reveals that younger clusters fit the 
observed line fluxes only when the upper mass cutoff' is smaller. 
For instance, a cluster younger than about 4 Myr can only re- 
produce the observed line ratios if its Mup is lower than 5OM0. 
An even younger cluster of about 1 Myr could only fit the ob- 
servations if its Mup is lower than 30Mq. 

The large spectral range covered by ISO included line 
fluxes such as [Arn] 7.0, [Nein] 15.5 and [S iii] 18.7 //m. These 
lines fluxes (listed in Tabled are upper limits to the true fluxes 
emitted by the nebula because of the large ISO aperture. They 
are compared to the model predictions in Fig.|3 As it was al- 
ready pointed out in Sect. 13.21 the comparison of the models 
with the observed line fluxes reveal that most of the emission 
from the high-excitation [Nem] Une (the ionisation potential 
of Ne^^ is 41 eV) comes from C2. However, other sources in- 
cluded in the ISO aperture must contribute to the low-excitation 
[Arn] and [Sm] line fluxes (Ar^ and S^^ have ionisation po- 
tentials of 15.8 and 23.3 eV) since the models predict much 
lower fluxes. Indeed, these fluxes must be taken as upper lim- 
its. Considering that practically all the [Ne m] line flux mea- 
sured by ISO is emitted by the embedded SSC C2, we can con- 
firm that the upper limit to the elemental abundance of neon 
estimated in Sect. l3.4l (0.4xrNe/H1iTi') is close to the true neon 
abundance of the H ii region. Since neon is a primary product 
of the stellar synthesis of elements and its abundance closely 



follows that of oxygen, we can affirm that the H ii region ex- 
cited by NGC 5253 C2 is metal-poor with a metallicity close to 
that measured for the galaxy. 

The effect of varying the ionisation parameter U (i.e. Rm 
and/or e) and other parameters such as the slope of the IMF, the 
nebular and stellar metal content, the internal dust content and 
the density structure, as well as the effect of a matter bounded 
geometry are discussed below. 

4.3.1 . The ionisation factor 

We have adopted 7?out = 0.8 pc, Ri„ = 0.6 pc and e = 1 for 
the models presented here. As mentioned in Sect. 14.21 possible 
solutions for the geometry of the nebula include lower values 
of Ri„ and e. Such solutions will produce an increase of the 
ionisation factor, as was shown in Fig. |2l and of the general 
ionisation level of the nebula. Hence, the models presented in 
Figs.|3]and|4]provide an upper limit to the upper mass cutoff for 
a given starburst age in the case of the nebula having a smaller 
inner radius (or a larger ionisation parameter). 

Another aspect that can be analysed with respect to the ion- 
isation factor is that of how it can be lower. A decrease of Qo or 
an increase of the size of the H n region will indeed lower the 
parameter U. Assuming that Qo is well constrained from the 
observed Bro- line flux, other geometries must satisfy the con- 
dition given by Eq. An upper limit to the size of the nebula 
is given by the 11.7 //m image of Goriianet al. (2001), where 
the source is seeing as an unresolved object with a diameter of 
0'.'58 (~ 9 pc). We wiU take then Rom - 4.5 pc as the largest 
possible outer radius of the nebula. Using this outer radius, the 
condition given by Eq.0requires an electron density larger than 
2830 cm"^. The lowest value for the ionisation parameter is 
obtained considering the largest possible values for the inner 
radius and electron density (see Eq. |3j- The electron density 
cannot be larger than that derived from the 7 mm observations 
(5 X 10^ cm"^) for a much smaller and compact Hii region. A 
nebula with Ri„ ~ Rout - 4.5 pc and n^-S x lO'* cm"^ will 
have an ionisation parameter equal to log U ~ -2.2. The filling 
factor of a nebular with this geometry is constrained by Eq. 
e(Rin^Rom)^O.l. 

We have run a set of photoionisation models using this ge- 
ometry («e=5 X 10"^ cm"-', Rin = 4.5 pc and e = 0.1) and stellar 
SEDs with Mup = lOOM© and Z = 0.008. These models are 
plotted by dashed-dotted-dotted lines in Figs. |3l and 0] They 
reproduce well the observed He i line flux and the [Ne m] mea- 
sured by ISO. Regarding the IR lines observed by TIMMI2, 
the [Ne ii] and [ Ar m] line fluxes are well reproduced by star 
clusters younger than 4 Myr, while the [S iv] line flux is under- 
predicted by a factor of about 4-15 in this age range. 

Stellar clusters with an upper mass cutoff of IOOMq can 
marginally reproduce all the observed line fluxes when the ra- 
dius of the Hn region is allowed to be as large as 4.5 pc. 
However, the Lyman continuum photon rate we estimate from 
the Bra line flux (~ 2 x 10^^ s ') is of the order of the value ob- 
tained from the 7 mm radio observations (~ 1 x 10^^ s"' when 
a distance of 3.25 Mpc is assumed) we use to constrain the size 
of the nebula. Hence, it does not seem plausible that the neb- 
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Fig. 3. Variation of the selected emission line ratios as a function of starburst age. The models reproduce the observed Bra line 
flux within 5%. The observed ratio is indicated by an horizontal line and its uncertainty by an error bar on the left side. Upper 
limits correspond to a mix of ISO and small aperture measurements (see TableQ. The legend of the figures is the following: lines 
in black correspond to models with Z=0.008; lines in light colour to models with Z=0.004; solid lines correspond to models with 
Mup - IOOMq; dashed lines to models with M^p = 50Mq and dotted lines to models with M^p = 30Mq. Dashed-dotted lines 
correspond to models with solar metallicity and Mup - 1OOM0. The nebular parameters of all these models are Rin = 0.6 pc. 
Rout = 0.8 pc, € - 1 and ne= 5 x lO'* cm"^. Dashed-dotted-dotted lines correspond to models with Mup = IOOMq, Z = 0.008 and 
the following nebular parameters: Rm ~ Rout - 4.5 pc, e = 0.1 and «£= 5 x 10"* cm""*. 




Fig. 4. Variation of the mid-infrared line ratios as a function of starburst age. The observed value is indicated by an horizontal 
line and its uncertainty by an error bar on the left side. Upper limits correspond to a mix of ISO and small aperture measurements 
(see TableQ. Models are plotted using the same code as in Fig.|3 



ula has a radius much larger than the value of 0.8 pc we have 
considered. 

4.3.2. The power law index of the IMF 

We have adopted a Salpeter IMF wit h a power la w index 
a - 2.35. However, for massive stars, IScald (Il986l) finds a 
steeper IMF with a = 2.7. Correcting for binaries can also lead 
to systematically larger exponents (a ^ 2.4 - 3; Kroupa 2002). 
A larger value (a ^ 3.0 + 0.1) is also suggested by a completely 
independent but indirect approach based on the FIR luminos- 
ity function of IRAS point sources with colours of ultracom- 

S iact H n regions and CS(2-1) detections in the Galactic disk 
Casassus et a l.'2000V 

We analyse here the effect of a steeper IMF on the pho- 
toionisation model results. In order to do so, we have run a set 
of photoionisation models with the reference parameters listed 
in Tableland stellar SEDs with M^p = 100 Mq, Z = 0.008 
and a = 3.0. The effect of a steeper IMF with a - 3.0 on 
the considered line ratios is not very significant. In the age in- 
terval 0-6 Myr, the predicted [ Ar iii] and [Ne ii] line fluxes in- 
crease by factors < 1.2 and < 2.1, respectively, while the [S iv] 
line flux decreases by a factor > 0.8. Hence, there is a soften- 
ing of the cluster SED because of the smaller number of very 
massive stars; however, the effect on the mid-IR lines is small. 
Concluding, an IMF with a steeper slope does not reconcile 
stellar clusters with M^p - lOOM© and age < 4 Myr with the 
observed line fluxes. 

The effect of a top-heavy IMF (i.e. with a smaller power 
law index) would be a hardening of the SED because of 
the increased number of very massive stars. Hence, the fidu- 
cial models presented in Figs. |3l and |3 (with a Salpeter in- 
dex, a - 2.35) provide an upper limit to the upper mass 
cutoff for a given starburst age in the case of an IMF with 
a < 2.35. There are only min or indications of c lusters with 
top-heavy IMF's. For instance, Figer et all (Il999al) find an in- 
dex a ~ 1.6 for the massive Arches and Quintuplet clus- 
ters situated near the Galactic centre. There also indications 
of a top-heavy IMF in the massive M82-F cluster located in 



the starburst galax y M82. The mass-tq-hght ratio derived to- 
wards this cluster jSmith & Gallaghe3l200l h is significantly 
smaller that the ratio expected from a normal IMF, indicating a 
deficit of low-mass stars and a top-heavy IMF. However, other 
well-studied and resolved massive clusters such as 30 Doradus 
JSelmanetal. 1999) and NGC3603 ( Eisenhauer et al. 192^ 
have an IMF power law index close to the Salpeter value. 

4.3.3. Metal enrichment 

The effect of a possible metal enrichment of the parental molec- 
ular cloud is discussed here. We have considered an H ii region 
with the same geometrical parameters as in the fiducial models 
(see Tabled but considering solar stellar (and nebular) metal- 
licity. The predictions of the photoionisation models are plotted 
in Figs.|3land|3by a dashed-dotted line. The general effect of 
increasing the metallicity is softening the SED and lowering 
the general ionisation level of the nebular gas. A young stellar 
cluster with Mup = IQOMq and solar metallicity can reproduce 
the observed line fluxes for ages < 4 Myr 

However, the upper limits for the elemental abundances es- 
timated in Sect. 13.41 (especially that for Ne/H, which must be 
close to the true neon abundance of the Hii region, see the 
discussion above) do not give evidences of a possible metal 
enrichment. This test is anyway useful to show how wrong as- 
sumptions on the metallicity can lead to wrong conclusions on 
the cluster properties. 

4.3.4. Internal dust 

The effect of adding internal dust in the nebula is to increase 
the absorption of ionising photons and to modify the shape of 
the ionising spectrum. The new CLOUDY version 96.00-beta 5 
presents a very detailed treatment of the grain physics, where 
the heating and cooling of the gas by grain photoionisation- 
recombination is determined consistently and the grain size dis- 
tribution is resolved. Resolving the size distribution can lead to 
significant changes in the emitted spectrum at IR wavelengths. 
CLOUDY allows the use of different grain size distribution 
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functions. We will consider here two of these distributions: the 
ISM type, which reproduces the observed ratio of total to selec- 
tive extinction (i.e. Ry - 3.1), and the Orion size distribution, 
which is deficient in small particles and reproduces the larger 
Ry observed in Orion. 

We have run several photoionisation models using a stel- 
lar SED with Mup = 1OOM0, 1 Myr and Z = 0.008, and the 
nebular parameters of the fiducial models (see Table |3}- We 
have considered internal dust following the ISM and Orion 
size distributions with an abundance 0.1 times their standard 
abundance. This is basically a scale factor used to multiply the 
stored grain opacities. For both types of dust, we have fixed go 
such that we reproduce the Bra line flux within 5% (we note 
that the flux of Bra have only been corrected for foreground 
extinction). 

A first result of these dusty models is that, in order to re- 
produce the observed Bra flux, we get a larger Qq with respect 
to the value obtained for the non-dusty models. In particular, 
go increases to about 4.0 x 10^^ s"' and 6.6 x 10^^ s"' for the 
models with ISM and Orion-like dust grains, respectively. With 
respect to the mid-IR line fluxes, the models give the following 
ratios: log([Neii]/Bra')=-2.7 and -2.3, log([Ariii]/Bra')=-1.3 
and -1.0, log([Siv]/Bra)= 0.5 and 0.7, log([S iv]/[Neii])=3.2 
and 3.0, and log([Ariii]/[Neii])=1.4 and 1.3. Hence, we see 
that the general ionisation level increases with respect to the 
non-dusty models, increasing even more the difference with the 
observed line fluxes. This is because the absorption cross sec- 
tion per H nucleon decreases for energies above 13.6 eV and 
consequently, the presence of dust favours the ionisation of He 
and other ions with respect to H. 

As a conclusion, the presence of internal dust does not rec- 
oncile the stellar clusters with Mup = IOOMq and age < 4 Myr 
with the observed line fluxes. The non-dusty models presented 
in Figs.|3land0]provide an upper limit to the upper mass cutoff 
for a given starburst age in the case of the presence of internal 
dust. 

4.3.5. Density profile 

We have considered so far photoionisation models with a con- 
stant electron density equal to 5 x 10"^ cm"^. We will analyse 
now the effect of an electron density profile which is allowed to 
vary as a power law of the form We = "o x (^/^in) where 
is the electron density at the illuminated face of the nebula and 
R is the nebular radius. We have run a set of models with differ- 
ent values of «o and y using a stellar SED with Mup=lOOM0, 
1 Myr and Z - 0.008. Predictions of these models are shown in 
TableE] where they are compared with the fiducial model (#1) 
and the observations. 

The first test, which comprises models #2 and #3, analy- 
ses the effect of a density gradient with an inner density equal 
to 5 X 10"* cm"^. A density gradient produces a slight increase 
of the [S iv]/[Ne ii] line ratio and a thickening of the ionised 
shell when compared with the fiducial model #1. The discrep- 
ancy with the observations is even higher. This increase in 
[S iv]/[Ne ii] is principally due to the fact that the [S iv] 10.5 ;um 
fine has a lower critical density (3.7 x 10"* cm"^) than the [Ne ii] 



12.8 fim line (6.1x10^ cm"-'). Hence, a density gradient favours 
the enhancement of the [S iv] line with respect to the [Ne ii] line 
(and also with respect to the [Ar iii] line, with a critical density 
equal to 2.1 x 10^ cm"^). 

The second test (models #4 and #5) studies the effect of 
decreasing the inner density, As noted in Sect. 14.21 an Hii 
region with Qo - 2 x 10^^ s"' and /?out = 0.8 pc does not 
allow values of the electron density lower than 3.7 x lO'* cm""*. 
As it is shown by Eq. the smallest H ii region with Qo - 
2 X 10^2 s-i and ^ 5 x 10^ cm would have an outer radius 
of ~ 3 pc (with Rin - 0.45 pc), about 4 times larger than what 
is observed at radio wavelengths. In terms of the line fluxes, 
these models predict higher fluxes and a general increase of 
the excitation of the gas (as measured by [S iv]/[Neii]). The 
discrepancy with the observations is, again, higher than that of 
the fiducial model. 

The third and final test (models #6 and #7) analyses the 
effect of increasing the inner density. We study two cases: 
no = 5 X 10^ and 10^ cm . We only consider constant den- 
sities. At these high densities, the thickness of the shell is very 
small and density gradients will not produce any effect on the 
predictions. These models predict an important decrease of the 
[S iv]/[Ne ii] ratio. This effect goes in the right direction to rec- 
oncile the models with Mup^lOOMo and the observations. The 
model with «□ = 5 x 10^ cm"^ even reproduces rather well the 
[S iv]/[Ne ii] line ratio; however, it fails to reproduce the indi- 
vidual fluxes of [Arm], [Siv] and [Neii]. The increase of the 
density up to 10^ cm"-' comes along with an important weak- 
ening of the [Siv] line flux. This weakening is also seen in 
the [Ar iii] line. On the contrary, the flux of the [Ne ii] slightly 
increases. This is because, at this high density, collisional de- 
excitation of the upper levels of the [Ar iii] and [S iv] lines is 
important, while the [Ne n] line, with a higher critical density, 
is much less affected. Hence, an increase of the density, even 
though leads to a lower ionisation level of the nebula, fails to 
reproduce the observed [Ar iii] and [S iv] line fluxes, while the 
[Ne ii] line flux, although higher, is still underpredicted. 

In conclusion, variations in the density and/or density pro- 
file do not reconcile the models with Mup^lOOM© and 1 Myr 
with the observations. The models presented in Figs. |3] and |4] 
provide an upper limit to the upper mass cutoff for a given 
starburst age in the case of power-law density profiles with in- 
ner densities lower than that assumed in the fiducial models 
(5 X 10"* cm"-'). We have also shown that even though an in- 
crease of the density leads to lower values of the [S iv]/[Ne n] 
ratio, such models largely underpredict the [Arm], [Siv] and 
[Neil] line fluxes. 

4.3.6. Matter bounded geometry 

The nebular models presented here have an ionisation (or radia- 
tion) bounded geometry, where the outer limit of the H ii region 
is defined by a hydrogen front. In this geometry, the H ii region 
is optically thick to the hydrogen ionising continuum and has 
absorbed nearly all of it. H ii regions can be also matter (or den- 
sity) bounded if the outer limit is set by the edge of the cloud. 
In this case, the H ii region is optically thin to the incident con- 
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Table 4. Comparison between predictions of models with different density profiles (see Sect. 14. 3. 5t and the observations. These 
models have been run using a stellar SED with Mup=100Mo, 1 Myr and Z - 0.008. They reproduce the observed Bra line flux 
within 5%. 



Model 
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log [Ariii]/Bry 


log [S ivl/Bry 


log [Ne iiJ/Bry 


log [Siv]/[Neii] 


^in (pc) 


AL* (pc) 


#1 


5 X 10-* 





-0.73 


0.68 


-2.11 


2.80 


0.6 


0.19 


#2 


5x lO'* 


1 


-0.76 


0.75 


-2.17 


2.92 


0.6 


0.26 


#3 


5 X lO'* 


2 


-0.80 


0.84 


-2.23 


3.08 


0.6 


0.41 


#4 


5 X 10^ 





-0.59 


1.14 


-2.00 


3.15 


0.45 


2.48 


#5 


5 X 10 '' 


1 


-0.40 


1.24 


-1.72 


2.95 


0.45 


41.3 


#6 


5 X 10' 





-0.52 


-0.38 


-1.35 


0.97 


0.8 


1.6 X 10-3 


#7 


10' 





-1.50 


-2.86 


-1.25 


-1.61 


0.8 


4.8 X 10^*^ 


Obs 






0.08 + 0.18 


0.68 + 0.15 


-0.52 ±0.14 


1.23 + 0.15 







* Thickness of the shell. 



tinuum. The effect of a matter bounded geometry will be then 
to cut off the external parts of the nebula, i.e., the lower excita- 
tion regions, producing a nebula with a higher ionisation level. 
Consequently, if the nebular were matter-bounded, the upper 
mass cutoff of the IMF would need to be revised downward in 
comparison to the models discussed above (Figs.|31and|3. 

4.4. On the [On] emission in NGC 5253 

The high-excitation [O iv] line at 25.9 jjm has been observed in 
active galactic n uclei and in some starburst galaxies including 
NGC 5253 (e.g. lGenzel eTal.. 1998: Lutz et al. 1998) . The ion- 
isation edge for the creation of O iv (54.9 eV) is just beyond 
the Hen edge where normal stars emit little flux and hence, 
the origin of such high-excitation line in starbursts has been 
unclear. Different excitation mechanisms such as weak AGNs, 
super-hot stars (e.g. WR stars), planetary nebulae and ionising 
shocks related to the starburst activity have been proposed by 
iLutz et aP (Il998l) . These authors favour ionising shocks as the 
most likely explanation for the presence of [O iv] emission in 
starburst galaxies. Schaerer & Stasinska ( 1999a) , on the other 
hand, propose that the emission of [Oiv] 25.9 yum observed 
with ISO in the two dwarf galaxies NGC 5253 and IIZw40 
is due to the presence of hot WR stars. Independently of the 
origin of the [Oiv] emission, the mere presence of such high 
energy source(s) in NGC 5253 calls for an examination of their 
possible effects on the photoionisation models presented above 
and on the robustness of the results. We shall briefly address 
this now. 

First we assume that the [Oiv] 25.9 fim emission ob- 
served in the ISO-SWS aperture originates from the cluster 
C2. Figure |5] shows the observed and predicted [Oiv]/Bra' 
line ratio from the starburst models discussed above (Figs. |3 
and 13. Manifestly, the predicted strength of the [Oiv] line is 
much lower than observed. In the best case (the model with 
Mup=lOOM0 and Z - 0.004), [Oiv]/Bra' reaches few percent 
(~ 3%) between ~ 3 to 5 Myr when WR stars are pre sent. 
In contrast to the models of Sch aerer & St asinska ( 1999a), the 
discrepancy between the observed and predicted [O iv] emis- 



sion found here is due to 1) softer stellar atmos phere mod- 
els (mostly WR spectra including line blanketing: 'Smith et alJ 
.2002 ), and 2) a reliable constraint on the ionisation parameter 
Can the models be reconciled with the observed [O iv ] line 
flux? 

In principle this could be possible by increasing the ion- 
isation parameter U - e.g. by adopting a smaller inner ra- 
dius and/or a density gradient and/or internal dust and/or a 
top-heavy IMF (see Sect. 14. 3> . However, as discussed previ- 
ously, such an increase of U would cause a larger discrepancy 
between the models with Mup=lOOM0 and the other mid-IR 
lines. In particular, the [Arm] line flux, which is barely re- 
produced by the fiducial models with Mup^lOOM© would be 
largely underpredicted. 

A more promising way to circumvent the shortfall of the 
observed [Oiv] line flux is by considering that the adopted 
spectral energy distribution of the cluster close to and beyond 
the He ii edge is not well defined. This could be justified for 
several reasons. First, X-emission in early typ e stars can alter 
their i onising spectra at h igh energ ies (see e.g.|Pa uldrach e t alJ 
EooitlSchaerer fc de KoteT..1997..MacFarlane et al.J994 and 
discussion therein). Second, other sources of high energy pho- 
tons as binaries or shock s responsible of soft X -ray emission 
could be present in C2 JSummers et alJE'004h . In this situa- 
tion, the ionising spectrum could be increased at high energy 
such as to enhance the [Oiv] emission without altering the 
ionisation structure (and hence the emission line strengths) of 
the other observed ions whose ionisation potentials are lower 
(IP ~ 30^0 eV). Indeed, after several ad-hoc modifications 
of the library SEDs beyond ~ 45 eV, we have verified that a 
spectral energy distribution with log22/2() ~ -1.8^ can re- 
produce both the observed [O iv] line flux and the other con- 
straints. This corresponds to a harder SED than typically ex- 
pected for starbursts (cf. Smith & Gallagher 2001). For in- 
stance, the hardest Starburst99 model used above (a Z = 
0.004 burst at 4 Myr) has log 62/60 - -2.7, whereas in 

^ Here Qo and Q2 stand for the ionising flux capable to ionise H and 
He^ respectively. 
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other models one typically has log ^2/ Go ~ -5 to -6. In 
fact, if the X-ray lum inosity determined for region #19 by 
ISummers et al ] (l2004 is attributed to C2, one has Lx/i'boi ~ 
10 considerably higher than the "normal" X- ray luminos- 
ity of early type stars (Chlebowski et alJll989l) . Invoking a 
harder SED than predicted by standard synthesis models would 
therefore seem justified. Finally the postulated hardness of 

log QiIQq 1-8 would correspond to an optical line ratio of 

I(HeIIA4686)/I(H/3) ~ 3%, a hardness not unu sual in metal- 
poor H 11 regions and Wolf-Rayet galaxies (cf. iGuseya et alJ 
l2000l). 
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Fig. 5. Variation of [Oiv]/Bra' as a function of starburst age. 
The observed ratio is indicated by an horizontal line and its 
uncertainty by an error bar on the left side. Models are plotted 
using the same code as in Fig.|3] 



From the above we see that is possible to construct mod- 
els able to reconcile the observed ISO [O iv] 25.9 ;um emission 
with the other observables for the cluster C2 without altering 
our conclusions on the age, IMF etc. Furthermore, the required 
harder ionising spectrum does not seem physically unplausi- 
ble. On the other hand, it is not excluded that part or all of the 
[O iv] line flux is not related to C2 and comes from another ob- 
ject included in the large ISO-SWS aperture. In this case, our 
observations and modelling require, for consistency, that this 
source shows a rather faint [Arm] 9.0 fim emission and a [S iv] 
10.5 fim comparable or less than C2. Mid-infrared observations 
with Spitzer should be able to test these predictions. 

5. Discussion 

5.1. Importance of high spatial resolution observations 

We have shown the importance of high spatial resolution ob- 
servations of extra-galactic star clusters. Only such observa- 
tions provide accurate measurements of the Une fluxes emitted 



by the cluster itself and avoid contamination by other clusters, 
the local ISM, etc. For instance, in the case of NGC 5253 C2, 
while most of the high-excitation [S iv] and [Ne iii] line fluxes 
measured by the large ISO aperture is emitted by the H 11 region 
associated to the cluster, other sources largely contribute to the 
low-excitation [Arii], [Siii] and [Neii] line fluxes. For exam- 
ple, the remaining presumably "diffuse" [Ne n] emission in the 
large ISO aperture amounts to ~ 80% of the total emission. 

This has implications not only in the interpretation of line 
fluxes in term of the properties (age, IMF, etc.) of the em- 
bedded cluster (e.g. Crowther et al. 1999; Thornley et al. 2000; 
iRigby & R ieke 2004), but also on the modelling of the SED 
(e.g. IVanzi & SauvagQi2004l) . where contamination by back- 
ground radiation from the underlying galaxy and other sources 
might lead to erroneous conclusions. 

High spatial resolution observations also provide impor- 
tant constraints on the nebular geometry (size, electron den- 
sity, etc.). The interpretation of nebular lines in terms of prop- 
erties of the ionising cluster depends on photoionisation models 
which have as free parameters the SED of the cluster (Myp and 
age), the local abundance and the ionisation parameter U. In 
particular, the parameter U, which depends on the geometry, is 
generally strongly underconstrained, leading to wrong conclu- 
sions. Along these Unes, we have shown how two SEDs with 
diff'erent M^p can lead to similar conclusions when two differ- 
ent geometries are considered (see Sect. 14. 3.11 . The same holds 
when the local abundance is not properly taken into account (cf. 
Sect. l43H . 

We have thus demonstrated that the spatial scale of 
the observations greatly determine the mid-IR appearance 
of NGC 5253. Similar aperture effects affecting both the 
line and PAH emission fluxes have been observed by 
ISiebenmorgen et aP (l2004l) in luminous IR galaxies. However, 
it is likely that "global" mid-IR (ISO or even the smaller Spitzer 
aperture) diagnostics give to a large extent correct answers con- 
cerning AGN/starburst diagnostics. This is supported by the 
comparison with optical diagnostics, which mostly agree with 
the mid-IR ones (c f. .Genzeletal. 1998; Laurent et al. 2000.; 
IPeeters et alJl2004l) . However, more detailed and subtle diag- 
nostics (IMF, age, etc.) are probably more prone to aperture 
effects. 



5.2. Constraints on the age and upper mass cutoff 

It has been claimed that the IMF upper mass cutoff is not 
well defined, but is gi ven by the stochastic nature of the IMF 
(cf. |Elmegreen"l997'). However, NGC 5253 C2 is massive 
enough (see Tab le |5) to be populated with stars up to Mup (cf. 
IWeidner & Kroupal2004 . 

For NGC 5253 C2, strong constraints on the geometry (i.e. 
Lyman photon rate, density, filling factor and inner/outer radius 
of the nebula) based on high spatial observations in the near-IR 
(HST and Keck), mid-IR (TIMMI2) and radio (VLA) lead to an 
ionisation parameter log U > -0.5 dex. This restriction on the 
parameter U allows us to conclude that a cluster with an upper 
mass cutoff of IOOMq can only reproduce the observed line 
fluxes if its age is about 5-6 Myr, above the upper Umit set by 
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the non presence of supernovae (~ 4 Myr). A cluster younger 
than about 4 Myr can fit the observations if its M^p < 50 Mq. 
A young cluster of 1 Myr would only be consistent with an 
even lower upper mass cutoff (around 30 Mq). This result is 
practically independent on the power law index of the IMF (cf. 
Sect. I4.3.2I I. the internal dust content (cf. Sect. I4.3.4I I and the 
density profile (cf. Sect. 14.3.51 . 

The photoionisation models allow higher values of the up- 
per mass cutoff for ages < 4 Myr only in the case of an H ii 
region much larger than what the radio continuum images im- 
ply (cf. Sect. 14. 3.11 or when a solar metallicity, for which there 
is no indication, is considered for the cluster (cf. Sect. l4.3.3t . 

In short, we are left schematically with two possible solu- 
tions due to an age degeneracy: 1) a young (< 4 Myr) cluster 
with a clear evidence for a "non-standard", low, upper mass 
cutoff of the IMF, or 2j a cluster of ~ 5-6 Myr with a "stan- 
dard" high upper mass cutoff. The existing age constraints have 
been described above (Sect. lTTl . Let us now discuss these two 
cases and their implications. 

5.2.1 . A young cluster deficient in massive stars? 

A young age ($ 3-4- Myr) would agree with the lack of super- 
novae signatures in C2, but it would imply a "non-standard" 
IMF with a low upper mass cutoff {M^p 5 50Mq). Indeed, 
so far there is a general consensus from various studies of in- 
dividual stars in clusters that stars of at least ~ 100 Mq are 
formed and that most likely M^p ~ 120-150Mo quite inde- 
pendently of metallicity and other "environmental" conditions 
Je.g. Massev & Hu nter 19 98, Figer et al. 1 999a, Kroupa 20 02, 
Weidner & Kroup al 12004 and reviews in lOiknore & Howelll 
1998. and iCorbelli et al .''2004^. M^p > 50Mq has also been 
found from integrated light studies of (giant) extra-galactic H ii 
regions, H ii g alaxies and alike o bjects at sub-sola r metallicities 
(cf. reviews of lLeithereill998l and*Schaere j2003l) . Even at high 
metallicity, recent direct s tudies of t he stellar content indicate a 
high/normal value of M^p jPindao et all2002h . Claims for a low 
upper mass cutoff of the IMF have been made repeatedly (e.g. 
poldader et al. 1997; Bresolin et al. 1999; Coziol et al. 2001; 
Rigbv & Riekell2004l) . However, these claims are essentially 
based on less direct diagnostics involving nebular line obser- 
vations, large apertures and/or objects with complex structures 
which render their interpretation more difficult. Here we find, 
for the first time, an indication for a "non-standard", low upper 
mass cutoff of the IMF for an individual massive cluster. If true, 
this implies that the IMF is different and not fully populated up 
to the largest masses in some regions/objects. What is responsi- 
ble for this deviation remains, however, unknown as apparently 
similar clusters with a normal IMF exist (cf. list in Table|5jl. Is it 
related to the cluster mass, to its strong gravitational potential, 
or dependent on other properties? Do dynamical processes and 
interactions play a role in establishing the IMF in such a case? 
A better understanding of the IMF and its origin are needed 
before these questions can be answered. 

What are the implications from a massive star deficiency? 
This obviously depends on the physical quantity of interest and 
on how general such a case is. Studies of additional hidden 



SSCs will obviously be of interest to examine if our result is 
also found for other objects. For example, we simply note that 
the bolometric luminosity emitted by a massive young (^ 4- 
5 Myr) cluster with Mup = 30 Mq is ~ 1 mag smaller per unit 
stellar mass than for a normal Mup = 100 M©. In other words 
the contribution of such a hidden cluster to the dust heating and 
IR flux is reduced by the same amount. 



5.2.2. An "old" hidden compact cluster? 

The solution of an "old" age of ~ 5-6 Myr raises several ques- 
tion. First, how is it possible to "contain" and hide such a 
compact cluster so long? With a much higher electron den- 
sity (~ 5 X 10^ cm"^) than what it is commonly found in 
typical giant extra-galactic Hii regions (~ 10^ - lO^* cm""*), 
NGC 5253 C2 should be overpressured compared to the sur- 
rounding ionised gas. At such high pressure, the region would 
expand rapidly to reestablish pressure equilibrium. Following 
IWood & Churchwelll lll989bll . we find that with an age of 5- 
6 Myr and under no other additional external pressure, C2 
would have expanded until a size of 22-24 pc in diameter. 
However, this size is more than 10 times what it is observed at 
radio wavelengths (1.6 by 0.6 pc, cf. Sect. 14. 1> . The most sim- 
ple explanation seems that the cluster mass is large enough to 
allow its own gravitational co nfin ement (sel f-gravitv), as sug- 
gested bv iTurner et all (l2003l) and lTan & McKee (2004) . 

While several other compact hid den clusters have 
been known for several year s ( Kobulnickv & JohnsonI Il999t 
'Neff^ & Ulvestedl l2000t iBeck et aL .200 2) they are generally 
thought to be very young with ages («c 1-4 Myr) mainly based 
on compactness arguments. In some cases also, compact ther- 
mal radio sources seem to show optical counterparts providing 
age estimates of ~ 0-4 Myr ( Whitmore & Zhang 2002). Here, 
however, our analysis indicates for the first time that such hid- 
den clusters could live longer. The "average" fraction of the 
lifetime massive stars spend buried in their parent molecular 
clouds can also be estimated by comparing the number of "hid- 
den" O stars to the total observed in the galaxy. The hidden 
cluster C2 is powered by ~ 2 x 10''^ Lyman continuum photons 
per second, which implies about 2000 07V equivalent stars. 
Over a 30" aperture, the total (extinction-corrected) Bry line 
flux is ~ 5 X 10"^" W cm"^, which is equivalent to 5 x 10^^ s"' 
and 5000 07 V equivalent stars. Based on these estimates, O 
stars in NGC 5253 might spend approximately 40% of their 
lives in the embedded phase. This percentage is much larger 
than the 1 0%-20% estimated for Galac tic ultracompact H ii re- 
gions (cf . 'Wood & Churchwell" 1989a') and the 15% estimated 
for young ultradense H ii regions in the young starburst galaxy 
He 2-10 bv lKobulnickv & JohnsonI ( 1 1999t) . 

If this scenario is true, it remains to be understood what de- 
termines the duration of the "hidden" phase of SSCs allowing 
for significant variations of the latter. Is it the cluster mass, con- 
centration, its gravitational potential, external effects or others? 
The only clue we have so far is that the cluster C2 analysed here 
is among the most massive and compact ones (see Table |5}- 
Future detailed studies of young massive clusters should be 
able to able to shed more hght into these questions. 
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On the other hand, the "old" solution of ~ 5-6 Myr is some- 
what puzzling given the apparent absence of supernovae in this 
region of NGC 5253 (cf. l4.U . Indeed, given the high luminos- 
ity/Lyman continuum production of C2 (see Table |5} one ex- 
pects a supernova ra te of at least 100 SN /Myr after approx- 
imately 3 Myr (e.g. ICervino et al1l2000l) or in other words, 
at least a total of several hundred past SN explosions at an 
age of ~ 5-6 Myr from stars with initial masses M ^ 35-40 
Mq. However, it seems difficult to understand how such a large 
number of supernovae could be "hidden" without giving raise 
to non-ther mal radio emission especiall y in such a high density 
region (cf. iChevalier & FranssoDl200lh . 

5.3. Comparison with previous mid-IR constraints on 
the IMF 

In the case of NGC 5253 we are facing an apparent "contra- 
diction". The ISO observations show that NGC 5253 is, with 
II Zw 40, the object with the highest excitation within the ISO 
extra-galactic sample (as measured by [Neiii]/[Neii]). For this 
reason this object is considered as "unproblematic", i.e. com- 
patible with a normal Salpeter IMF extending to at least ~ 40- 
60 Mq ( Rigby & Rie ke 2004) in contrast to the majority of the 
starbursts observed by ISO. In addition, our observations now 
indicate an even higher excitation for the SSC C2. Why then do 
we here face this difficulty of reconciling the observations with 
a "normal", high value of M^p? The answer to this apparent 
contradiction is simple: the difficulty arises from the fact that 
a proper geometrical constraint is used here for the first time 
implying a fairly large ionisation parameter. Had we neglected 
the information on the spatial extent (or better said the com- 
pactness) of the emission associated with C2, we would easily 
be able to reconcile the observations with a standard value of 
Mup for any burst age. This case illustrates the potentially se- 
vere limitations inherent in the traditional interpretations of fine 
structure line ratios when the source distribution and geometry 
are unknown. However, from the single case study presented 
here it is not possible to generalise the conclusions. 

5.4. Comparison with other massive clusters 

A stellar cluster with M^p- 100 M© and an age of 6 Myr would 
have a total cluster mass Mdustei ~ 34x 10'' Mq and a total num- 
ber of stars A^stars =9x10^ for a Salpeter IMF slope extending 
down to 1 Mq. a stellar cluster with M^p- 50 Mq and an age of 
4 Myr would have Mduster ~ 9 X 10^ Mq and A^.tars ~ 2 x 10^ 
Finally, a stellar cluster with Mup = 30 Mq and 1 Myr would 
have Mciuster = 4 x 10^ M0 and N,,.,,, = 1 x 10^ 

Table |5l compares NGC 5253 C2 and other massive clus- 
ters in mass, size, density, age, luminosity and Lyman contin- 
uum flux. NGC 5253 C2 is most comparable in mass to MGG- 
11 and [W99]2, but it is younger. In terms of mass density, 
NGC 5253 C2 is comparable to the Galactic Centre and Arches 
clusters, but it has a higher luminosity and produces more ion- 
ising photons. 

Regarding the upper mass cutoff, studies of R136, the ion- 
ising cluster of 30 Doradus, find that the most probable upper 



mass limit in the Large Magellanic Cloud is about 130Mf ri (cf. 
Sel man et alj [l999). Towards the Arches cluster, Figer et alJ 
j 1999 a) find that there more than 10 stars with masses larger 
than 12OM0. If true, the case of a low upper mass cutoff in 
NGC 5253 C2 would greatly differ with the high values found 
in these massive clusters. 

5.5. Host of an intermediate mass black hole? 

Recent dynamical simulations bvl Portegies Zwart et alJ (l2004h 
show that massive stars in young dense clusters can sink rapidly 
to the centre of the cluster to participate in a runaway collision 
and produce a star of 800-3000 M0 which ultimately colla pses 
to a b lack hole of intermediate mass (see also Giirkan e t alJ 
I2OO4I and references therein). 

The conditions for the formation of a black hole inside a 
dense cluster depend on two distinct parameters, the concen- 
tration parameter (c) and the dynamical friction timescale (fdf). 
The parameter c is defined as c = logiRJRc), where Rt and 
Rc represent the tidal radius and the core radius, respectively; 
fdf is the time taken for a star of a certain mass to sink to the 
cluster centre from a circular or bit at in itial distance R sj> Rc . 
According to the simulations bv lPorteg ies Zwart et al.] ( l2004l) . 
the requirement for the formation of such a black hole is that 
the cluster is born with c > 2 and that fjf < 4 Myr The 
4 Myr upper limit for t^f is such that massive stars can eas- 
ily reach the centre of the cluster before exploding as super- 
novae. These authors' simulations indicate that in the cluster 
MGG-11, located in M82 and associated with an ultralumi- 
nous X-ray source, a runaway star could form on the necessary 
timescale. In the case of NGC 5253 C2, we find fdf ~ 2.5 Myr, 
calculated fo r a star of 50 Mn star t ing at a radius of 0.8 pc. 
As stated by Porte gies Zwart et al ] (|2a04'), less massive stars 
undergo weaker friction and thus must start at smaller radii in 
order to reach the cluster centre on a similar timescale. Hence, 
an intermediate mass black hole might form through runaway 
collisions in the core of NGC 5253 C2. 

Recent Chandra and XMM -Newton observations of 
NGC 5253 dSummers etalJl2004 have found a X-ray point 
source at the peak of the radio emission associated with C2 
(their source #19). However, this X-ray source is very soft, with 
very few counts above 1.0 keV and virtually none above 2.0 
keV. The authors suggest that this X-ray emission is associ- 
ated with the emission from superbubles surrounding the inner 
cluster rather than with an individual harder object such as a 
black hole. Besides, they find that the X-ray emission coinci- 
dent with C2 seems to be extended compared to the Chandra 
PSF, reinforcing the idea that the emission is associated with 
the cluster Hence, although it seems that X-ray observations 
show no evidence for the presence of a black hole within C2, it 
is a possibility that deserves certain consideration. 

6. Conclusions 

We have presented the A^-band (8-13 fim) spectrum of the hid- 
den compact radio super-star cluster in NGC 5253, C2, ob- 
tained with TIMMI2 on the ESO 3.6 m telescope. The spec- 
trum is characterised by a rising continuum due to warm dust, 
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Table 5. Properties of massive clusters. Given are the total cluster mass, radius, mass density in stars, age, luminosity and rate of 
Lyman photons. 





lUg IVlf IVlQ 


Kadius (pc) 


iogp(,MQpc ; 
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NGC 5253 C2 


< 6.5" 


0.8 


< 6.2 
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9.0 


52.3 


1, 2, 3, 4 


Quintuplet 
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1.0 


3.2 


3-6 
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50.9 
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0.19 
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1-2 


8.0 


51.0 


5 


Gal. Center 


4.0* 


0.23 


5.6 


3-7 


7.3 


50.5 


5 


NGC 3603 


3.7* 


0.23 


5.0 
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7.3 


51.1 


5 


R136 


4.5* 


1.6 


3.3 


< 1 -2 


> 7.6 


51.9 


5 


MGG-11^ 


5.5^ 


1.2 


4.7 


7-12 






6 


[W99]2* 


6.3'' 


4.5 
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" If due to virialised motions, the Bry line width measured bv iTurner et alJ j22Q3) indicates a cluster mass of the order of 10* M© for a 
Salpeter IMF slope and a lower mass cutoff of ~ 1 Mq or less. * Total cluster mass in all stars extrapolated down to a lower mass cutoff 
of 1 Mq, assuming a Salpeter IMF slope and an upper mass cutoff of 120 M ^. ^ Kinematic mass. ^ Star cluster in M82. ^ Star clust er in 
NGC 4038/4039. R EF ERENCES: (1) This wo rk; (2) Turner & Beck (2004); r3) IVanzi & SauvageH2004h ; f4l lAlonso-Herrero et alJliool : (5) 
lFigeretalJil999d) ; (6) lMcCradv etai] <2003h ; ('7) iMengel et ah (.2002.1 . 



a silicate absorption and a strong [S iv] line at 10.5 jim. Weaker 
lines of [Arm] at 9.0 yum and [Ne ii] at 12.8 yum are also present. 
The continuum can be modeled by an optically thick emission 
from hot (T^ - 253 + 1 K) dust emission extinguished by a 
cold foreground dust screen and a silicate absorption feature 
with Asii = 0.73 + 0.05 mag. 

It has been demonstrated that the spatial scale of the obser- 
vations greatly determine the mid-IR appearance of NGC 5253. 
In particular, most of the high-excitation [Arm], [Siv] and 
[Ne III] line fluxes measured by the large ISO aperture are emit- 
ted by C2, while other sources largely contribute to the low- 
excitation [Aril], [Siii] and [Nen] line fluxes. This has impor- 
tant implications in the interpretation of line fluxes in terms of 
the properties (age, IMF etc.) of the embedded cluster. For in- 
stance, the present data indicate an even higher excitation for 
C2 compared to the observations through the large ISO/SWS 
aperture. 

We have computed sets of nebular models with the pho- 
toionisation code CLOUDY using the evolutionary synthesis 
code Starburst99 to model the integrated properties of the stel- 
lar cluster. It has been shown that the interpretation of the neb- 
ular lines in terms of properties of the ionising cluster depends 
largely on the local abundance and the ionisation parameter 
U. The detailed dependence of the mid-IR lines on other pa- 
rameters such as the cluster age, upper mass cutoff and power 
law index of the IMF, as well as the presence of internal dust 
and the density structure is largely discussed. In the case of the 
SSC C2, high spatial resolution observations at different wave- 
lengths - near-IR (HST and Keck), mid-IR (TIMMI2) and ra- 
dio (VLA) - have allowed us to strongly constrain the geom- 
etry of the region (i.e. Lyman photon rate, density, filling fac- 
tor and inner/outer radius), leading to an ionisation parameter 
log U > -0.5 dex. This constraint on U lead to two possible 
solutions for the age and upper mass cutoff of C2 when com- 
paring the observed line fluxes with the predicted ones: 7) a 
young (< 4 Myr) cluster with a "non-standard" IMF with a low 
upper mass cutoff Mup< 50 Mq, and 2) a cluster of ~ 5 - 6 Myr 



with a standard high upper mass cutoff (Mup~ 100 Mq). The 
photoionisation models allow higher values of Mup for ages 
< 4 Myr only in the case of an H ii region much larger than 
what the radio continuum imply or when a solar metallicity, 
for which there is no indication, is considered for the cluster A 
young age (< 4 Myr) would agree with the lack of supernovae 
signatures in C2 and in case of being confirmed, would be the 
first indication for a "non-standard", low upper mass cutoff of 
the IMF for an individual massive cluster An older age of ~ 
5-6 Myr would imply that it is possible to "contain" and hide 
such a compact cluster for a longer time that what it is generally 
thought. Arguments in favour and against these two scenarios 
are presented. Based on the available information, we are not 
able to favour one solution over the other. A better understand- 
ing of the signature of supernovae from massive star progeni- 
tors in compact and dense environments, as well as radio ob- 
servations at even higher spatial resolution, would permit to 
exclude or not the ages above 4 Myr Besides, mid-IR observa- 
tions at even higher spatial resolution (which can be performed 
e.g. in the near future with VISIR on the VLT), would allow 
determining the exact size of C2 at these wavelengths and re- 
solving its ionisation structure. In particular, resolving the ion- 
sation structure of C2 would be useful to verify if the lines we 
are using for our photoionisation analysis come from the same 
region or if there is any contamination by the close presence 
of the optical stellar cluster CI. Observations of massive stellar 
clusters with similar characteristics will be as well essential to 
see if the same result is found. 

We have addressed as well the origin of the [Oiv] 25.9 jum 
emission measured by ISO assuming that all of it originates 
from the SSC C2. The strength of the [Oiv] line predicted 
by the photoionisation models is much lower than observed. 
However, it has been shown that it is possible to construct 
models able to reconcile the observed ISO [Oiv] emission 
with the other mid-IR lines by modifying the adopted SED 
of the cluster close to and beyond the Hen edge. A SED 
with log22/2i ~ -1-8 can reproduce all the observables. 
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This hardness would correspond to an optical line ratio of 
I(HeIIA4686)/I(H/3) ~ 3%, not unusual in metal-poor Hii re- 
gions and Wolf-Rayet galaxies. 

Finally, based on dynamical simulations, we consider the 
possibility that an intermediate mass black hole might form 
through runaway collisions in the core of C2. However, recent 
X-ray observations show that the X-ray source associated with 
C2 is very soft and it is most probably associated with the emis- 
sion from superbubles surrounding the inner cluster. 
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Appendix A: The 1 8 jum flux of NGC 5253 C2 

Besides the A^-band spectrum, we also obtain an image of 
NGC 5253 C2 in the Ql band (50% band cuts are 17.35 and 
18.15 fim). The object appears unresolved, with a FWHM of 
2'.'4 + O'.'l. The flux was calibrated using the standard star 
HD 123139. We obtain a total flux of 4.53 + 0.25 Jy. This 
val ue is ~ 1.6 times l arger than the 18.7 fim flux reported 
bv lGoriian et alJ (1200 Ik 2.9 + 0.3 Jy. However, according to 
Gorjian (private communication), their low 18.7 jim flux could 
be due to a too small chopping throw, resulting in an incorrect 
zero-level, or to a complex airmass coiTection. 

References 

Aflamandola, L. J., Tielens, A. G. G. M., & Barker, J. R. 1989, 
ApJS, 71,733 

AUende Prieto, C, Lambert, D. L., & Asplund, M. 2001, ApJ 
Lett., 556, L63 

Alonso-Herrero, A., Takagi, T, Baker, A. J., et al. 2004, ApJ, 
in print 

Audouze, J. 1987, in lAU Symp. 124: Observational 

Cosmology, Vol. 124, 89-115 
Beck, S. C, Turner, J. L., Ho, P. T P, Lacy, J. H., & Kelly, 

D. M. 1996, ApJ, 457,610 
Beck, S. C, Turner, J. L., Langland-Shula, L. E., et al. 2002, 

AJ, 124, 2516 

Boker, T, Krabbe, A., & Storey, J. W. V. 1998, ApJ Lett., 498, 
L115 

BresoHn, F, Kennicutt, R. C, & Garnett, D. R. 1999, ApJ, 510, 
104 

Campbefl, A., Terlevich, R., & Melnick, J. 1986, MNRAS, 223, 
811 

Casassus, S., Bronfman, L., May, J., & Nyman, L.-A. 2000, 
A&A, 358,514 

Cervino, M., Knodlseder, J., Schaerer, D., von Ballmoos, P., & 
Meynet, G. 2000, A&A, 363, 970 



Charbonnel, C, Meynet, G., Maeder, A., Schafler, G., & 

Schaerer, D. 1993, A&AS, 101, 415 
ChevaHer, R. A. & Fransson, C. 2001, ApJ Lett., 558, L27 
Chlebowski, T, Harnden, F R., & Sciortino, S. 1989, ApJ, 341, 

427 

Cohen, M., Walker, R. G., Carter, B., et al. 1999, ApJ, 117, 
1864 

CorbelU, E., Palla, F, & Zinnecker, H. 2004, IMF@50, The 
Initial Mass Function 50 years later (Astrophysics and Space 
Science, in press) 
Coziol, R., Doyon, R., & Demers, S. 2001, MNRAS, 325, 1081 
Crowther, P A., Beck, S. C, WilUs, A. J., et al. 1999, MNRAS, 
304, 654 

Eisenhauer, F, Quirrenbach, A., Zinnecker, H., & Genzel, R. 

1998, ApJ, 498, 278 
Elmegreen, B. G. 1997, ApJ, 486, 944 

Ferland, G. J., Korista, K. T, Verner, D. A., et al. 1998, PASP, 
110, 761 

Figer, D. F, Kim, S. S., Morris, M., et al. 1999a, ApJ, 525, 750 
Figer, D. F, McLean, I. S., & Morris, M. 1999b, ApJ, 514, 202 
Freedman, W. L., Madore, B. F, Gibson, B. K., et al. 2001, 
ApJ, 553, 47 

Gurkan, M. A., Freitag, M., & Rasio, F A. 2004, ApJ, 604, 632 
Garay, G., Rodriguez, L. F, Moran, J. M., & Churchwell, E. 

1993, ApJ, 418, 368 
Genzel, R. & Cesarsky, C. J. 2000, ARA&A, 38, 761 
Genzel, R., Lutz, D., Sturm, E., et al. 1998, ApJ, 498, 579 
Gilmore, G. & Howell, D. 1998, The Stellar Initial Mass 

Function, 38th Herstmonceux Conference (ASP Conference 

Series, Vol. 142) 
Goldader, J. D., Joseph, R. D., Doyon, R., & Sanders, D. B. 

1997, ApJ, 474, 104 
Gorjian, V, Turner, J. L., & Beck, S. C. 2001, ApJ Lett., 554, 

L29 

Guseva, N. G., Izotov, Y. I., & Thuan, T X. 2000, ApJ, 531, 
776 

Hiflier, D. J. & Miller, D. L. 1998, ApJ, 496, 407 
Home, K. 1986, PASP, 98, 609 

Hummer, D. G. & Storey, P J. 1987, MNRAS, 224, 801 
Kobulnicky, H. A. & Johnson, K. E. 1999, ApJ, 527, 154 
Kobulnicky, H. A., Skillman, E. D., Roy, J., Walsh, J. R., & 

Rosa, M. R. 1997, ApJ, 477, 679 
Kroupa, P 2002, Science, 295, 82 

Kurucz, R. 1993, ATLAS9 Stellar Atmosphere Programs and 
2 km/s grid. Kurucz CD-ROM No. 13. Cambridge, Mass.: 
Smithsonian Astrophysical Observatory, 1993., 13 

Laurent, O., Mirabel, I. F, Charmandaris, V, et al. 2000, A&A, 
359, 887 

Leitherer, C. 1998, in ASP Conf. Ser. 142: The Stellar Initial 
Mass Function (38th Herstmonceux Conference), 61 

Leitherer, C, Schaerer, D., Goldader, J. D., et al. 1999, ApJS, 
123,3 

Lumsden, S. L., Puxley, P J., & Doherty, R. M. 1994, MNRAS, 
268, 821 

Lutz, D., Kunze, D., Spoon, H. W. W., & Thornley, M. D. 1998, 

A&A, 333, L75 
MacFarlane, J. J., Cohen, D. H., & Wang, P 1994, ApJ, 437, 

351 



Martm-Hemandez et al.: Mid-IR spectroscopy of the embedded SSC in NGC 5253 



19 



Madden, S. C. 2002, Ap&SS, 281, 247 

Martm-Hemandez, N. L., van der Hulst, J. M., & Tielens, 

A. G. G. M. 2003, A&A, 407, 957 
Martin-Hernandez, N. L., Peelers, E., Morisset, C, et al. 2002a, 

A&A, 381, 606 

Martm-Hernandez, N. L., Vermeij, R., Tielens, A. G. G. M., 

van der Hulst, J. M., & Peeters, E. 2002b, A&A, 389, 286 
Massey, R & Hunter, D. A. 1998, ApJ, 493, 180 
Mathis, J. S. 1990, ARA&A, 28, 37 

McCrady, N., Gilbert, A. M., & Graham, J. R. 2003, ApJ, 596, 
240 

Mengel, S., Lehnert, M. D., Thatte, N., & Genzel, R. 2002, 
A&A, 383, 137 

Meurer, G. R., Heckman, T. M., Leitherer, C, et al. 1995, AJ, 

110, 2665 

Meynet, G., Maeder, A., Schaller, G., Schaerer, D., & 

Charbonnel, C. 1994, A&AS, 103, 97 
Neff, S. G. & Ulvestad, J. S. 2000, AJ, 120, 670 
Okamoto, Y. K., Kataza, H., Yamashita, T., Miyata, T., & 

Onaka, T. 2001, ApJ, 553, 254 
Okamoto, Y. K., Kataza, H., Yamashita, T., et al. 2003, ApJ, 

584, 368 

Osterbrock, D. E. 1989, Astrophysics of gaseous nebu- 
lae and active galactic nuclei (Research supported by 
the University of California, John Simon Guggenheim 
Memorial Foundation, University of Minnesota, et al. Mill 
Valley, CA, University Science Books) 

Pagel, B. E. J. 1992, in lAU Symp. 149: The Stellar Populations 
of Galaxies, Vol. 149, 133 

Panagia, N. & Walmsley, C. M. 1978, A&A, 70, 41 1 

Pauldrach, A. W. A., Hoffmann, T. L., & Lennon, M. 2001, 
A&A, 375, 161 

Peeters, E., Martm-Hemandez, N. L., Damour, P., et al. 2002, 
A&A, 381,571 

Peeters, E., Spoon, H. W. W., & Tielens, A. G. G. M. 2004, 

ApJ, in print; astro-ph/0406183 
Pindao, M., Schaerer, D., Gonzalez Delgado, R. M., & 

Stasiriska, G. 2002, A&A, 394, 443 
Portegies Zwart, S. P., Baumgardt, H., Hut, P., Makino, J., & 

McMillan, S. L. W. 2004, Nature, 428, 724 
Rigby, J. R. & Rieke, G. H. 2004, ApJ, 606, 237 
Roche, P R & Aitken, D. K. 1984, MNRAS, 208, 481 
Rubin, R. H., Simpson, J. P., Erickson, E. P., & Haas, M. R. 

1988, ApJ, 327, 377 
Salpeter, E. E. 1955, ApJ, 121, 161 

Scalo, J. M. 1986, Pundamentals of Cosmic Physics, 11,1 

Schaerer, D. 2003, in in lAU Symposium 212, 642 

Schaerer, D., Charbonnel, C, Meynet, G., Maeder, A., & 

Schaller, G. 1993a, A&AS, 102, 339 
Schaerer, D. & de Koter, A. 1997, A&A, 322, 598 
Schaerer, D., Meynet, G., Maeder, A., & SchaUer, G. 1993b, 

A&AS, 98, 523 
Schaerer, D. & Stasiriska, G. 1999a, A&A, 345, L17 
Schaerer, D. & Stasiriska, G. 1999b, in ESA SP-427: The 

Universe as Seen by ISO, 751 
Schaller, G., Schaerer, D., Meynet, G., & Maeder, A. 1992, 

A&AS, 96, 269 
Sehnan, P., Melnick, J., Bosch, G., & Terlevich, R. 1999, A&A, 



347, 532 

Siebenmorgen, R., Kriigel, E., & Spoon, H. W. W. 2004, A&A, 

414, 123 

Smith, L. J. & Gallagher, J. S. 2001, MNRAS, 326, 1027 
Smith, L. J., Norris, R. R R, & Crowther, R A. 2002, MNRAS, 

337, 1309 

Soifer, B. T., Neugebauer, G., Matthews, K., Egami, E., & 

Weinberger, A. J. 2002, AJ, 124, 2980 
Soifer, B. T., Neugebauer, G., Matthews, K., et al. 2001, AJ, 

122, 1213 

Summers, L. K., Stevens, I. R., Strickland, D. K., & Heckman, 

T. M. 2004, MNRAS, 351, 1 
Tan, J. C. & McKee, C. P. 2004, in Cancun Workshop on 

Rormation and Evolution of Young Massive Clusters, eds. 

H. J. G. L. M. Earners, A. Nota and L. Smith, see astro- 

ph/0403498 

Thomley, M. D., Schreiber, N. M. R, Lutz, D., et al. 2000, ApJ, 

539, 641 

Thuan, T. X., Sauvage, M., & Madden, S. 1999, ApJ, 516, 783 
Tumer, J. L. & Beck, S. C. 2004, ApJ Lett., 602, L85 
Turner, J. L., Beck, S. C, Crosthwaite, L. P., et al. 2003, Nature, 

423, 621 

Tumer, J. L., Beck, S. C, & Ho, R T. R 2000, ApJ Lett., 532, 
L109 

Turner, J. L., Ho, P T. P, & Beck, S. C. 1998, AJ, 1 16, 1212 
Vacca, W. D., Johnson, K. E., & Conti, R S. 2002, AJ, 123, 772 
Vanzi, L. & Sauvage, M. 2004, A&A, 415, 509 
Verma, A., Lutz, D., Sturm, E., et al. 2003, A&A, 403, 829 
Walsh, J. R. & Roy, J. 1989, MNRAS, 239, 297 
Weidner, C. & Kroupa, P 2004, MNRAS, 348, 187 
Weingartner, J. C. & Draine, B. T. 2001, ApJ, 548, 296 
Whitmore, B. C. & Zhang, Q. 2002, AJ, 124, 1418 
Wood, D. O. S. & Churchwell, E. 1989a, ApJ, 340, 265 
Wood, D. O. S. & Churchwell, E. 1989b, ApJS, 69, 831 



